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FOREWORD 


This  report  is  divided  into  two  completely  separate  and  distinct  parts. 
In  Part  I,  Physical  Optics  Inverse  Scattering t  a  physical  optics  approximation 
to  the  monostatic  electromagnetic  inverse  scattering  problem  is  developed  for 
stationary  perfectly  conducting  scatterers.  In  Part  II,  Exact  Inverse 
Scattering ,  a  rigorous  and  exact  solution  to  the  general  inverse  problem  is 
developed,  which  is  applicable  to  the  monostatic-bistatic  electromagnetic 
inverse  scattering  problem  for  moving  scatterers  of  arbitrary  conductivity. 
This  division  into  two  distinct  parts  is  necessitated  by  the  totally 
different  approaches  taken  to  the  problem;  although  the  results  of  Part  I 
are  reaffirmed  by  the  incidental  first-order  approximations  developed  in 
Part  II. 


ABSTRACT 


PART  I:  A  three-dimensional  electromagnetic  Inverse  Scattering  Identity, 
based  on  the  Physical  Optics  approximation,  is  developed  for  the  monostatic 
scattered  far  field  cross  section  of  perfect  conductors.  Uniqueness  of  this 
inverse  scattering  identity  is  proven.  This  identity  requires  complete 
scattering  information  for  all  frequencies  and  aspect  angles.  A  non-singular 
integral  equation  is  developed  for  the  arbitrary  case  of  incomplete  frequency 
and/or  aspect  angle  scattering  information.  A  general  closed  form  solution 
to  this  integral  equation  is  developed,  which  yields  the  shape  of  a  scatterer 
from  such  incomplete  information.  A  specific  practical  radar  solution  is 
presented.  The  resolution  of  this  solution  is  developed;  yielding  short- 
pulse  target  resolution  radar  system  parameter  equations.  Results  of  the 
three-dimensional  numerical  reconstruction  of  a  sphere  and  a  cylinder  from  a 
variety  of  aspect  angle  and  frequency  band  limited  cross  section  data  are 
presented.  The  merits  of  this  solution  over  the  conventional  synthetic 
aperture  radar  imaging  technique  are  discussed. 

PART  II:  The  general  inverse  scattering  and  radiation  problem 
associated  with  the  three-dimensional  inhomogeneous  scalar  field  Helmholtz 
wave  equation  is  formulated  as  a  Fredholm  integro-differential  equation  of 
the  second  kind.  The  far-field  inverse  integro-differential  equation  is 
solved  in  closed  form  with  the  aid  of  a  single  resolvent  integral  operator, 
which  can  be  readily  evaluated  numerically  with  the  aid  of  the  fast  Fourier 
transform  algorithm.  The  inverse  integro-differential  equation  and  its 
solution  are  then  generalized  to  the  reduced  vector  wave  equation  resulting 
from  Maxwell’s  equations.  A  formal  statement  of  the  inverse  problem  is 
presented.  It  is  shown  that  the  first  order  Neumann  series  solution  of  the 
inverse  integro-differential  equation  as  well  as  the  first  order  term  of 
its  exact  solution  represent  the  physical  optics  approximation  and  the 
equations  governing  synthetic  microwave  holography.  The  (analogous)  four¬ 
dimensional  inverse  integro-differential  equation  and  its  closed  form 
solution,  applicable  to  Doppler-contaminated  fields  radiated  by  moving 
scatterers,  is  developed. 
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PHYSICAL  OPTICS  INVERSE  SCATTERING 


SECTION  I 


INTRODUCTION 


The  Direct  Scattering  problem,  whether  electromagnetic,  acoustic, 
particle,  or  quantum  mechanical,  is  defined  as  the  problem  of  predicting 
the  scattered  quantities,  given  the  incident  quantities,  the  relevant 
description  of  the  scatterer,  and  the  appropriate  laws  governing  the  inter¬ 
action.  The  back -scattering  and  forward-scattering  problems  are  the  special 
monostatic  case  and  the  special  bistatic  case  of  scattering  in  the  direction 
of  incidence  respectively. 

The  Inverse  Scattering  problem  is  defined  as  the  problem  of  determining 
the  relevant  quantities  describing  the  scatterer,  given  the  incident  and 
scattered  quantities,  and  the  appropriate  laws  governing  the  interaction. 

The  Uniqueness  and  Well-behavedness  problems  must  also  be  taken  as  an  integral 
part  of  the  inverse  scattering  problem.  The  Well-behavedness  problem  is  the 
problem  of  determining  the  degree  of  smoothness  and  continuity  with  which  the 
so-called  output  data  varies  with  respect  to  the  so-called  input  data.  The 
uniqueness  problem  must  also  deal  with  the  question  of  incomplete  input  data. 

The  Electromagnetic  Inverse  Scattering  problem  is  thus  defined  as  the 
problem  of  determining  the  size,  shape  and  electromagnetic  properties 
distributions  (conductivity,  susceptibility,  and  permeability  distributions) 
of  a  scatterer,  given  the  incident  and  scattered  electromagnetic  fields,  and 
the  electromagnetic  field  equations  (Maxwell's  equations  and  the  appropriate 
wave  equations;  the  constitutive  equations  and  their  coefficients  are  taken 
as  part  of  the  electromagnetic  properties  of  the  scatterer);  and  the  determi¬ 
nation  of  whether  this  problem  is  uniquely  solvable  for  incomplete  input 
data,  i.e.,  the  various  permutations  of  incomplete  bistatic  aspect  angles, 
incomplete  monostatic  aspect  angles,  incomplete  frequency,  monochromatic  data 


i 
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only,  incomplete  polarization  matrix,  amplitude  (power  cross-section)  data 
only,  and  scattered  far  field  data  only.  Also  of  interest  is  the  inverse 
scattering  problem  for  which  some  a  priori  information  about  the  scatterer 
exists;  e.g.,  that  the  conductivity  of  the  scatterer  is  infinite  (perfect 
scatterer),  that  the  geometry  of  the  scatterer  is  of  axial  symmetry,  etc. 

For  a  physically  more  sophisticated  and  mathematically  more  rigorous 
definition  of  the  general  inverse  problem  (including  the  inverse  radiation 
problem),  the  reader  is  referred  to  part  II,  sections  IX  and  X.5  of  this 
report . 


The  subject  matter  of  this  part  I  is  restricted  to  the  monostatic  far 
field  special  case;  a  restriction  being  characteristically  inherent  to  the 
radar  application.  For  mathematical  reasons  that  will  become  evident  to  the 
reader,  the  further  special  case  of  a  priori  knowledge  of  the  scatterer  being 
a  perfect  conductor  is  treated.  For  both  of  these  cases,  special  attention 
is  given  to  the  further  special  case  of  the  finite  three-dimensional  aperture; 
i.e.,  the  case  for  which  monostatic  scattering  data  is  available  for  a 
limited  and  incomplete  frequency  and  aspect  angle  domain.  From  a  practical 
short  pulse  target  resolution  radar  point  of  view,  this  special  case  is  of 
fundamental  and  primary  interest. 

i 

The  approach  taken  and  the  results  obtained  can  be  summarized  as 
follows.  In  Section  II,  the  physical  optics  approximation  is  taken  as  the 
basis  for  the  direct  scattering  theory,  and  its  validity  for,  and  consistency 
with,  short  pulse  target  resolution  radar  concepts  is  discussed  briefly. 

In  Section  III,  with  the  aid  of  this  approximation,  a  basic  Inverse  Scattering 
Identity  is  then  developed,  that  states  that  the  characteristic  function  (in 
three-dimensional  space)  of  a  scatterer  is  related  to  the  field  cross-section 
(in  three-dimensional  k-space)  by  a  three-dimensional  Fourier  Transform. 
Uniqueness  of  the  solution  for  finite  sized  scatterers  is  established.  In 
Section  IV,  a  general,  non-singular,  inverse  scattering  integral  equation  is 
developed,  solutions  to  which  permit  the  determination  of  an  appropriate 
maximum  of  information  about  a  scatterer  from  incomplete  scattering  data. 

The  use  of  the  Physical  Optics  approximation  is  further  justified.  In  Section 


V,  a  closed  form  solution  to  this  integral  equation,  valid  for  any  arbitrarily 
shaped  finite  k-space  aperture  (of  incomplete  frequency  and  aspect  angles 
scattering  data),  is  developed.  The  full  details  of  a  practical  frequency  and 
aspect  angles  band  limited  right  rectangular  quasi-conic  section  aperture  are 
then  developed  as  an  example.  In  Section  VI,  the  spatial  resolution  of  the 
solution  obtained  is  developed,  yielding  a  set  of  short  pulse  target  resolution 
radar  system  parameters.  In  Section  VII,  the  special  cases  of  two-  and  one¬ 
dimensional  inverse  scattering,  the  special  case  of  certain  given  simple 
geometrical  shapes  of  the  aperture,  and  the  special  case  of  a  priori  knowledge 
of  certain  symmetries  possessed  by  the  scatterer,  are  discussed.  The  special 
case  of  a  priori  knowledge  of  the  scatterer  being  a  surface  of  revolution 
about  some  axis  is  treated  in  detail.  In  section  VIII,  results  of  the  three- 
dimensional  numerical  reconstruction  of  a  sphere  from  a  variety  of  aspect  angle 
and  frequency  band  limited  cross  section  data,  obtained  numericaly  by  Mie's 
exact  solution  for  the  direct  scattering  by  a  sphere,  are  presented.  Similar 
results  for  a  cylinder  are  also  presented.  In  section  IX,  plans  for  future 
numerico-experimental  work  are  outlined,  and  the  extension  of  the  monostatic 
physical  optics  inverse  scattering  solution  for  stationary  perfectly  con¬ 
ducting  scatterers  to  a  generalized  monostatic-bistatic  inverse  scattering 
solution  for  moving  scatterers  of  arbitrary  (unknown)  conductivity  is  briefly 
discussed.  In  the  concluding  section  IX,  conclusions  are  drawn  about  the 
merits  of  the  presented  solutions  over  the  conventional  synthetic  aperture 
radar  imaging  technique. 
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SECTION  II 


THE  PHYSICAL  OPTICS  CROSS  SECTION 

s 

The  scattered  magnetic  field  H  in  terms  of  the  induced  (hy  the  incident 
field)  surface  current  density  K  on  the  surface  of  a  perfect  conductor  is  given 
by  [1] 


=  < )  V$xK  ds 
•  /* 


CD 


where  the  Green's  function  <J>  and  its  gradient  are  given  by 

i  kr 


*  = 


A-nr 


_ .  _  I  kr- 1  ^  , 
V4>  =  — psr  r  *  • 


(2) 

(3) 


The  Physical  Optics  Approximation  [2]  for  the  induced  surface  current 
density  in  terms  of  the  incident  magnetic  field  is 


K  = 


ZnxH',  on  the  "illuminated"  segment  of  s, 
0  ,  on  the  "shadow"  segment  of  s. 


(4) 


f 

Thus,  by  (1),  (3),  and  (4),  the  scattered  magnetic  far  field  W  ,  in  terms 
of  the  wave  number  propagation  vector  k  of  the  scattered  far  field,  is 


=  -2iksx  (j,  nxH^  ds  . 
k^*n<0 


(5) 
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If  the  incident  field  n  is  taken  as  a  plane  wave  of  the  form 


H  =  I  e 


ik^-x 


in  the  vicinity  of  the  scatterer,  where  k  is  the  wave  number  propagation 
vector  of  the  incident  field*  and  the  range  and  phase  normalized  (ir  the 
coordinate  system  in  which  the  scatterer  is  described)  scattered  far  field 


is  taken  as 


- - S  e 

/4Vr2 


Ik  *x 


then  for  the  case  of 


It8  =  -jt1,  =  It  , 


i.e.,  the  monostatic  case,  (5)  reduces  with  the  aid  of  (2),  (8),  and  the 
transvexsality  of  the  incident  field  to 


(.-**•«  k- 

k*n>0 


ds  I  . 


Consistent  with  the  conventional  definition  [3]  of  the  power  cross  section 
o  and  the  field  cross  section  p 


S  =  p  I 


in6!2 

o  =  PP*  =  4irr2  — —  , 

jff|2 


equation  (9)  yields  for  the  physical  optics  field  cross  section  the  well- 
known  expression  [4] 


e-2lk'x  k*ds  . 
k*n>0 


For  a  short-pulse  target  resolution  radar  system  to  be  effective,  its 
pulse  length  must  be  short  compared  to  the  target  size;  furthermore,  since 
the  fractional  bandwidth  of  such  a  pulse  is  limited  by  practical  considerations 
to  much  less  than  unity,  it  follows  that  the  largest  wave  length  in  the 
spectrum  of  the  transmitted  pulse  must  be  very  short  indeed  compared  to  the 
target  size.  The  physical  optics  approximation  (12)  is  thus  a  good  model 
for  the  direct  scattering  theory  for  such  a  short-pulse  target  resolution 
radar  system.  For  a  detailed  discussion  of  the  physical  meaning  and 
implications  of  the  physical  optics  approximation,  i.e.,  its  being  a  total 
first  order  local  scattering  theory,  consistent  with  short-pulse  radar 
concepts,  the  reader  is  referred  to  an  earlier  work  of  this  author  [5]  . 


SECTION  III 


THE  INVERSE  SCATTERING  IDENTITY 


Introducing  the  variable  k  defined  as 


k  s  2k 


(13) 


yields  for  (12) 


p(K) 


K*n>0 


x*ds  . 


(14) 


Thus 


*  /  \ 
p  (-K) 


K*n<0 


(15) 


and 


p(K)  +  p  *  ( -K ) 


'  i  -!K*X  „  .  '  f  -iK*X  „ 

-  I  e  K*ds  +  — -  I  e  k< 

v^iT  J  Sat!  J 

K-n>0  K*n<0 

1  -Ik*x  „ 

— —  9  e  <c*dS  . 

JAn  J 


ds 


(16) 


(17) 


SscsDiic  ?*os  bu*-*t 


v  Since  the  integTand  of  (17)  is  continuous  and  differentiable  on  s  and  in 
v  banded  by  s,  it  follows  by  Gauss'  theorem  that 


pOO  +  p  <-«)  = 


V-Oce"'***)  dv 


K 


i 


dv  . 


(18) 


(19) 


Introducing  the  quantity  F(k)  defined  by 


r(K)  =  V^ir 


pOO  +  p  (-«) 


,2 


(20) 


yields  for  (19) 


f  ~1K.) 

J6 


r(K)  =  |  e  "  dv  . 
v 


(21) 


Defining  the  characteristic  function  y(X)  of  the  scatterer  by 


y(x)  s 


1  ,  x  in  v  , 

0  ,  X  not  in  v  , 


(22) 


permits  the  reformulation  of  (21)  as  the  three-dimensional  Fourier  integral 


roe 


y(X)  d3x  . 


(23) 


If  the  volume  v  of  the  scatterer  is  finite,  then  by  (22) 


1 


|y (x) |  d3x  =  v 


<  oo 


(24) 


10 


it  thus  follows  from  (23)  that  for  finite  sized  scatterers  the  three-dimensional 
inverse  Fourier  transform  of  r C*)  exists  uniquely;  i.e. 


With  the  aid  of  (20}  »  (25)  can  be  reformulated  as  follows 


(25) 


y  (x) 


1 

4TT5/2 


P*(-K) 


(26) 


which,  after  replacing  K  t>y  in 


the  second  integral  of  (26)  yields 


y(x)  = 


(27) 


(28) 


(29) 


Both  inverse  scattering  identities  (25)  and  (29)  clearly  require 
complete  scattering  information;  namely,  knowledge  of  p(K)  over  all  K-space 
(i.e.,  all  frequencies  an«3  3I4.  aspect  angles). 
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SECTION  IV 


TOE  FINITE  APERTURE  INTEGRAL  EQUATION 

In  practice  p(K)  is  known  (.measurable)  only  for  an  incomplete  finite 
portion  of  the  complete  /c-space;  namely  a  fc-space  aperture  consisting  of  a 
limited  (finite l  frequency  band  and  a  limited  aspect  angles  band.  Kurthermoi 
(23)  is  valid  only  in  the  Physical  Optics  regime  (wave  length  short  compared 
to  the  overall  size  of  the  scatterer) ,  and  hence  (25)  must  either  include 
fictitious  (Physical  Optics  scattering  data  in  the  Rayleigh  regime,  which  is 
physically  not  realizable)  low  frequency  scattering  data  r(JO,  or  no  such 
data  at  all.  It  is  thus  to  the  problem  of  determining  what  can  be  deduced 
about  a  scatterer  (i.e.,  y(X))  from  such  limited  high  frequency  finite 
aperture  data  that  the  ensuing  sections  are  addressed. 

Let  A(k)  be  an  aperture  function  defined  as 

A (k)  =  C(k)  W OO  ,  (30) 

where  C(k)  is  a  characteristic  aperture  function  defined  as 


COO 


1  ,  for  k  for  which  J'(k)  is  known, 

0  ,  for  k  for  which  r(K)  is  unknown, 


(31) 


and  where  W(K )  is  any  appropriately  chosen  (in  general  non-zero)  aperture 
weighting  function ,  subject  to  the  conditions 


[ 


| A(K '  \  d3<  <  00 


(32) 
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Thus,  if  the  *c-space  volume  of  the  aperture  is  finite,  and/or  the 
aperture  weighting  function  is  appropriately  chosen,  then  the  three- 
dimensional  inverse  Fourier  transforms  of  the  aperture  and  characteristic 
aperture  functions  exist  uniquely,  i.e. 


a(x) 


1 

(2it)  3 


c(x)  = 


(2tt)  3 


A(k)  d3K 


C U)  d 3k  . 


(33) 


(34) 


Thus,  by  (25)  and  the  three-dimensional  convolution  theorem  for  three- 
dimensional  Fourier  transforms,  it  follows  that 


a(x)*y(x)  = 


r(K)  A («)  d3<  , 


which  by  (30)  and  (31)  reduces  to 


a(x)*y(X)  = 


r(/c)  W(K)  d3<  , 


(35) 


(36) 


where  r(K)  is  clearly  known  in  the  domain  of  integration  C:  i.e,,  the 
aperture.  The  right-hand  side  of  (36)  can  thus  be  taken  as  known;  say  the 

known  function  g(X);  i.e. 


g(x)  sT2^T3  eiK*X  r(k)  W(K)  d3<  ’  (37) 

'L 

The  three-dimensional  inverse  scattering  problem  for  a  finite  aperture 
thus  reduces  by  (36)  and  (37)  to  the  three-dimensional  non-singular  convolution 
integral  equation  (a  Fredholm  integral  equation  of  the  First  Kind) 
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a(X-X’)  yix')  d3x'  =  g(X) 


This  integral  equation  (38)  can  be  solved  numerically  by  a  variety 
of  existing  techniques  such  as  the  matrix  methods  of  Ritz-Galerkin  [6] ,  the 
associated  Least  Square  Best  Estimate  method  [7] ,  and  the  associated  moments 
method  of  Harrington  [8] ,  the  Eigen-function  expression  method  of 
Toraldo  Di  Francia  [9],  leading  to  so-called  super-resolution,  and  the 
k-space  method  of  this  author  [10],  which  also  leads  to  super  resolution. 
Several  closed  form  solutions  of  (38)  for  apertures  of  specific  geometry 
have  been  obtained  by  Lewis  [11] ;  an  alternate  closed  form  solution  of  (38) 
for  apertures  of  general  arbitrary  geometry  is  presented  next . 


SECTION  V 


A  SOLUTION  OF  THE  FINITE  APERTURE  INTEGRAL  EQUATION 


The  solution  of  (38)  for  y(X)  is  greatly  facilitated  by  the  special 
properties  of  y(X)  (i.e.,  a  priori  knowledge  that  y(X)  is  a  characteristic 
function  of  the  form  (22)  )  and  the  possible  judicious  choice  of  the  aperture 
function  F(* ) . 

Let  the  x3-axis  be  chosen  as  passing  through  the  (near)  center  of 
the  aperture  A  (see  fig.  1).  Next,  let  the  aperture  function  W(k)  be 
chosen  as 


WOO  =  i*3  (39) 

It  thus  follows  from  (30),  (32),  (33),  and  (34),  and  again  the  three- 
dimensional  convolution  theorem  for  three-dimensional  Fourier  transforms  that 

a(X)  =  (.:(X)*6fx1)6(X2)6' (x3)  ,  (40) 


where  c(x)  is  known,  i.e. 


c(x)  = 


(27 r)3 


d3< 


(41) 


It  thus  follows  from  (38)  that 

c(X)*6(x1)6(x2)6'(x3)*y(X)  =  g(X)  ,  (42) 

which  reduces  to 
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K3>*3 


FIGURE  1 

THE  APERTURE  AND  SURFACE  COORDINATE  SYSTEM 
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c(x)»  |^3  =  g(x)  . 


(43) 


Since  y*x)  is  a  characteristic  function  of  the  form  (22);  namely  and 
particularly,  that  for  fixed  values  of  xj  and  x2>  y(x3)  is  a  dual  step 
'unction  in  X3  of  unity  magnitude  and  steps  at  the  lower  and  upper  surfaces, 
f.ay  Z'(x1,x2)  and  Z(xi,x2)  respectively,  of  the  scatterer  (see  fig.  1);  it 
'•allows  that 


-  6(x3  -  Z' (xi ,x2) )  -  6(x2  -  Z(xj,x2))  ,  (44) 

which  by  (43)  yields 

c(x)*6(x3  -  Z' (xj ,x2) )  -  c(X)*6(x3  -  Z(x],x2))  =  g(X)  .  (4^) 

Examination  of  (20) ,  (26),  and  (37),  and  symmetry  and  physical 
considerations  (i.e.,  the  implications  of  (4)  and  (12)),  thus  yields  for 
the  upper  surface  Z(xi,x2)  only 


c(x)*6 (x3  -  Z(x1.,x2)) 


provided 


W(K)  =  W*(-K)  , 


(46) 


(47) 


which  is  assured  by  (39) . 


It  is  physically  reasonable  (by  the  implications  of  (4)  and  (12)  ) 

that  information  about  the  lower  surface  Z’Cxi,x2)  should  only  be  obtain- 

£ 

able  from  scattering  data  from  the  lower  image  aperture  A  (-*>.  It  is  now 
evident  that  the  introduction  of  the  image  aperture  served  the  sole  purpose 
of  a  mathematical  artifice  which  permitted  the  application  of  Gauss'  theorem 
to  (12),  yielding  (25)  laid  (29);  and  that  knowledge  of  scattering  data  in 
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this  image  aperture  is  not  needed. 


The  three-dimensional  convolution  on  the  left-hand  side  of  (46),  say 
X(X),  reduces  to 


X<X)  =  c(X)#6(x3  -  Z(xj,x2)) 


(48) 


00 

JJ 


c(xi-xi,x2-x^,x3-z(xl,xp)  dxj  dx£ 


(49) 


Thus 


(50) 


where  x<X)  is  a  resolution  density  function  which  is  a  measure  of  the  location 
of  the  upper  surface  Z(xj,x2).  That  x^  is  indeed  such  a  resolution  density 
function  can  best  be  visualized  by  considering  the  limiting  case  of  an 
infinite  aperture  function  A(k)  for  which  c(X)  =  6(x);  for  such  an  aperture, 
by  (49) 


XCX)  = 


00 

■I 


6<x  -x[) 


6<x  -x’> 
2  2 


6(x3-Z(xJ,x^))  dxj 


dx^ 


(51) 


=  6(X3  -  Z(xj,x2))  ;  (52) 

whereas  for  a  practical  realistic  aperture  of  finite  K-space  extent,  the 
spatial  extent  of  the  non-vanishing  portion  of  c(X),  and  hence  the  non¬ 
vanishing  portion  of  x(x),  is  still  small  compared  to  the  size  of  the 
scatterer  (see  Sect.  VI).  In  fact,  it  is  this  resolution  function  x<X) 
which  determines  the  resolution  of  the  solution  (50);  a  resolution  which  can 
only  be  exceeded  by  the  super-resolution  method  mentioned  earlier. 
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A  three-dimensional  density  plot  of  |x(X)|  thus  represents  the 
smeared  geometrical  image  of  the  surface  of  the  scatterer;  the  spatial  extent 
of  the  smearing  clearly  being  the  spatial  extent  of  c(X),  i.e.,  the 
resolution. 

A  Best  Estimate  of  can  alternatively  be  obtained  by  a  variety 

of  correlation  (between  (41)  and  (49))  methods,  employing  Fourier  transform 
theory. 


SECTION  VI 


TOE  RESOLUTION 


The  resolution  in  x- space  is  clearly  the  spatial  extent  of  the  nearly 
non-vanishing  extent  of  y(X).  It  thus  follows  from  (30)  through  (.38)  and  (49) 
that  the  resolution  in  any  one  direction  in  x-space  is  the  reciprocal  of  the 
K-spatial  extent  of  the  non-vanishing  portion  of  the  Characteristic  Aperture 
Function  COO  in  that  same  one  direction. 

The  finite  aperture  inverse  scattering  solution  (50)  can  clearly  be 
reformulated  in  a  variety  of  desired  practical  (radar)  spherical  coordinate 
systems.  For  the  particular  spherical  coordinate  system  shown  in  fig.  2, 
i.e. 


K 


(cos  £ 
sin  £ 
cos  £ 


sin  c 
cos  £ 


(53) 


d3<  =  k?  d<  cos  £  d£  d£ 


(54) 


the  inverse  scattering  solution  (50)  for  the  right  rectangular  quasi-conic 
section  aperture  shown  becomes 


X<x)  = 


4tt5/2 


<2  Co  £o 

nitcxo  sin  £  I 

J 

Kl  “Co  “£ o 


1  it  cos  £  (xj  sin  £  +  X2  cos  £) 


p(x>£>t) 


X  cos  £  d£  cos2£  d£  k  dtc  ,  (55) 


where 
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where  k0  and  Ate  are  the  carrier  and  bandwidth  of  the  transmitted  spectrum 
respectively. 

Examination  of  (41)  and  (49),  after  a  similar  reformulation  into  this 
same  spherical  coordinate  system,  readily  reveals  (after  small  angles 
approximations)  the  resolution  in  range  and  cross  ranges  to  be 


4X1  *  K.  «  : 

A5  =  2?0 

(58) 

a*2  *  K.  «  : 

AC  =  250 

(59) 

ax3  -  • 

(60) 

This  resolution  can  only  be  exceeded  by  the  earlier  mentioned 
techniques  of  super  resolution. 

Equations  (58),  (59),  and  (60)  are  thus  the  equations  for  the  parameters 
of  a  target  resolution  radar  system. 


SECTION  VII 


TOE  SPECIAL  CASES 


The  special  cases  of  the  one-  and  two-dimensional  inverse  scattering 
problem  (i.e.,  scattering  data  restricted  to  a  K-space  line  or  plane 
respectively,  obtained  from  a  three-dimensional  scatterer)  can  be  treated  by 
applying  the  methods  of  Lewis  [12]  or  of  this  author  [13]  to  (30)  et  seq. 

Namely,  by  choosing  the  characteristic  aperture  function  for  the  two-  and 
one- dimensional  special  cases  respectively  as 


COO 

=  C(ki,k2)  <5(k3) 

(61) 

coo 

=  C(kj  )  6  (k2)  6  (k3) 

(62) 

The  three-dimensional  convolution  integral  equation  (38)  reduces  respectively 
to  the  two-  and  one-dimensional  integral  equations 

a(xi,x2)#6(xi,x2>  =  g(xj,X2)  (63) 

a(xi)*a(xj)  =*  g(xi>  (64) 

where  6  and  a  are  the  thickness  distribution  function  in  the  x3-direction 
and  the  area  distribution  function  orthogonal  to  the  xi-direction  of  the 
scatterer  respectively,  and  g(xi,x2)  and  g(xj)  reduce  respectively  to  (see 
(37)  et  seq.) 


g(xi>><2)  = 


\*1+k2x2) 


r(K1,K2)  W(k1,k2)  d*i  ck2 


(65) 
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r(<i)  W(ki)  cIkj 


(66) 


g(xi>  =  2'!-feI,q><1 

JC 


The  non-singular  integral  equations  (63)  and  (64)  can  be  solved  for  3 
and  a  by  any  of  the  previously  discussed  means. 


The  further  special  case  of  a  priori  knowledge  of  the  scatterer  being 
a  surface  of  revolution  about  the  x2-axis  (see  fig.  2)  can  clearly  be  treated 
by  the  preceding  two-dimensional  formulation  (63)  and  (65)  by  recognizing 
that  the  profile  function  (generatrix  of  revolution)  of  the  scatterer  is 
{B(0,x2);  thereby  further,  simplifying  (63)  and  (65)  after  the  appropriate 
modifications.  A  more  direct  treatment  of  the  problem  of  the  surface  of 
revolution  will  be  presented  subsequently. 


The  special  cases  of  the  aperture  A(k)  being  of  certain  given 
geometrical  shapes  can  be  treated  by  applying  the  method  of  Lewis  [14]  to 
(30)  et  seq. 

The  special  cases  of  a  priori  knowledge  of  the  scatterer  possessing 
certain  geometrical  symmetries  can  be  treated  by  applying  the  methods  of  this 

author  [15]  to  (30)  et  seq.  As  an  illustrative  example,  the  case  of  the 
scatterer  known  to  be  a  surface  of  revolution  about  the  X2~axis  (see  fig.  2), 
is  presented  next. 

For  such  a  scatterer,  the  monostatic  cross-section  clearly  is 
independent  of  the  longitudinal  aspect  angle  s;  i.e. 

p(k,C,C>  =  p(ic,5)  .  (67) 

Furthermore,  the  profile  function  X3(x2)  of  such  a  surface  of  revolution 
(see  fig.  2)  is  given  by  the  function  describing  the  upper  surface  Z(xi,X2> 
at  the  plane  xj=o.  It  thus  follows  from  taking  the  limits  of  integration 
over  K  in  (50)  as  from  0  to  2n,  setting  xi=o  in  (50),  and,  with  the  aid  of 
the  integral  representation  of  the  Bessel  functions,  that 


1 


Ji(kx3  aos  C)  p(k>£)  cos2£  d£  k  die  , 


(68) 


X(X!,X3)  = 


2h3/2 


<2  So 

If- 

<1  “£o 


Ikxj  sin  ? 


which  for  small  angle  approximations  yields  the  doubly  truncated  two- 
dimensional  mixed  Hankel  (Fourier-Bessel) -Fourier  transform 


X<*i,x3) 


<2  So 

<1  “Co 


IkCxi 


p(k,£)  dC  k  die  . 


It  is  noteworthy  that,  as  expected  intuitively,  only  two-dimensional 
scattering  information  (in  k  and  £;  i.e.,  in  frequency  and  one  aspect  angle 
the  latitudinal  aspect  angle)  is  required  for  an  inverse  scattering 
solution  by  (68)  or  (69) . 


Solutions  are  also  obtainable  for  the  various  combinations  of  the 
special  cases  cited. 


(69) 


SECTION  VIII 


NUMERICAL  RESULTS 


The  solution  to  the  integral  equation  (38)  proposed  by  Lewis  [11]  was 
successfully  numerically  tested  for  a  sphere  by  this  author  in  1969.  This 
test  consisted  of  a  computer  implementation  of  a  special  case  version  of  this 
solution,  applicable  only  to  scatterers  about  which  only  a  priori  knowledge 
of  cylindrical  symmetry  exists.  This  test  essentially  confirms  the  correctness 
of  the  basic  inverse  scattering  identity  (25)  and  the  finite  aperture  integral 
equation  (38).  This  solution  was,  however,  not  pursued  further  because  of  its 
inherent  practical  limitations.  These  limitations  are  the  lack  of  generality 
of  the  required  k-space  aperture  (i.e.,  the  required  aperture  is  impractical 
for  physically  realizable  radar  systems;  which  is  not  the  case  with  this 
author’s  solution  50),  the  error  enhancement  introduced  by  the  process  of 
numerical  differentiation  of  noisy  scattering  data  (vis-a-vis  the  error 
reduction  resulting  from  the  process  of  integration  of  such  data  in  solution 
50) ,  and  the  unapplicability  of  the  Fast  Fourier  Transform  (FFT)  to  this 
solution  (which  is  essential  if  large  amounts  of  data  are  to  be  processed  in 
reasonable  time  by  existing  computers,  yielding  three-dimensional  high- 
resolution  descriptions  of  arbitrarily  shaped  scatterers  about  which  no  a  priori 
knowledge  of  special  geometry  exists) . 

Solution  (50)  was  computer  implemented  with  the  aid  of  the  FFT  for 
arbitrarily  shaped  apertures,  realizable  with  existing  radar  systems.  This 
computer  program  was  tested  with  the  exact  solution  of  Mie  for  scattering  by  a 
sphere,  with  a  variety  of  band  limited  aspect  angles  and  fractional  frequency 
band  widths  8  (see  fig.  3),  with  the  results  shown  in  fig.  5  through  8.  The 
basic  inverse  scattering  identity  was  also  tested  by  this  program;  with  results 
shown  in  fig.  4. 
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FIGURE  7 


RESOLUTION  DENSITY  FUNCTION  X  AT 


g'  =  48.36  % 
U  *  14.48° 
x3  »  o 


FIGURE  8 


x(*i)  x2 


RESOLUTION  DENSITY  FUNCTION  x  AT 


3'  -  9.00  % 

U  -  2.61° 

x3  *  o 


3. 


Only  X3=0  "slices"  of  the  reconstructed  scatterer  are  shown  in  fig.  4 
through  8;  however,  several  other  than  x3=0  "slices"  were  also  numerically 
reconstructed  correctly,  thus  confirming  the  three-dimensional  capabilities 
of  solution  (50) .  A  full  band-width  three-dimensional  display  of  a  recon¬ 
structed  sphere  and  cylinder  are  shown  in  fig.  9  and  10  respectively. 


FIGURE  9 

THREE-DIMENSIONAL  DISPLAY  OF  RECONSTRUCTED  SPHERE 


FIGURE  10 

THREE-DIMENSIONAL  DISPLAY  OF  RECONSTRUCTED  CYLINDER 


SECTION  IX 


PRESENT  STATUS  AND  FUTURE  WORK  HANS' 


The  fast.  Fourier  transform  computer  program  utilized  in  the  numerical 
work  summarized  in  sect.  8  was. an  in-core  FFT  program  limited  to  about 
16,000  data  points.  For  increased  resolution,  particularly  at  high  fre¬ 
quencies,  low  band  width,  and  small  angular  apertures,  it  is  necessary  to 
have  an  FFT  program  operative  on  about  (100) 3  =  106  data  points.  Since  data 
of  this  size  cannot  be  stored  in  core,  and  must  be  stored  externally  (e.g., 
drum  or  disc),  an  external- storage- FFT  program  is  currently  being  implemented 
under  a  follow-on  contract.  Under  this  follow-on  contract,  the  inverse 
scattering  solution  utilizing  such  an  external-storage  FFT  will  also  be 
tested  against  actual  measured  scattering  data  from  a  variety  of  modeled 
shapes  of  varying  complexity. 

The  physical  optics  inverse  scattering  solution  developed  in  this 
part  I  has  been  extended  from  the  presented  monostatic  solution  for 
stationary  perfectly  conducting  scatterers  to  a  generalized  monostatic- 
bistatic  inverse  scattering  solution  for  moving  scatterers  of  arbitrary 
(unknown)  conductivity.  This  generalized  physical  optics  solution  turns  out 
to  be  identical  with  the  first  order  solution  of  the  exact  inverse  scattering 
solution  presented^ in  part  II.  It  is  for  this  reason  that  this  generalization 
is  not  presented  here. 
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SECTION  X 


CONCLUSIONS 

The  solutions  developed  in  this  part  I  can  be  viewed  as  solutions  to 
the  frequency  band  limited  and  aspect  angles  limited  short  pulse  synthetic 
aperture  radar  imaging  (and  associated  data  processing)  problem.  The 
solutions  presented  are  based  on  rigorous  electromagnetic  scattering  and 
inverse  scattering  theory  applicable  to  spatially  distributed  scatterers, 
yielding  real  three-dimensional  geometrical  images;  vis-a-vis  the 
conventional  two-dimensional  radar  imaging  technique  which  is  based  on  the 
heuristic  approach  of  isolated  point  scatterers  (scattering  centers) 
concepts,  yielding  so-called  radar  images  (or  maps)  [16],  which  most  often 
do  not  resemble  the  real  geometrical  images  of  the  scatterer. 

The  merits  of  the  solutions  presented  in  this  paper  over  the 
conventional  radar  imaging  technique  deserve  the  following  further  exami¬ 
nation.  The  rather  unsatisfactory  results  of  the  conventional  radar  imaging 
technique  are  essentially  the  consequence  of  the  technique  being  based  on  a 
heuristic  approach  to  the  problem;  i.e.,  the  approach  consists  of  considering 
a  spatially  extended  target  as  a  fictitious  ensemble  of  identifiable, 
stationary,  non- interactive,  non-dispersive,  and  isotropic  point  scatterers. 
From  a  rigorous  electromagnetic  scattering  point  of  view,  a  spatially 
extended  scatterer  is  neither  an  ensemble  of  point  scatterers,  nor  are  these 
fictitious  point  scatterers  in  principle  identifiable,  stationary,  non¬ 
interactive,  non-dispersive,  and/or  isotropic  scatterers.  (The  point 
scatterers  are  not  always  identifiable  by  virtue  of  the  so-called  registration 
problem;  i.e.,  that  the  point  scatterers  can  be  continuously,  consistently, 
and  correctly  identified  for  various  aspect  angles.  The  point  scatterers 
are  not  always  stationary  due  to  the  fictitious  relocation  caused  by  changing 
aspect  angles.)  Furthermore,  this  technique  does  not  contain  basic  unique 
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existence  considerations,  and/or  rigorous  considerations  of  the  problem  of 
optimizing  the  results  for  incomplete  aspect  angles  and/or  frequency  infor¬ 
mation  availability.  The  attempts  to  convert  radar  images  to  geometrical 
images  have. thus  failed  for  precisely  these  reasons. 

The  solutions  of  this  part  I  alleviate  all  these  objections  to  the  so- 
called  radar  imaging  technique  by  the  rigorous  application  of  electromagnetic 
inverse  scattering  theory,  based  on  rigorous  direct  scattering  theory  (and 
not  based  on  the  heuristic  model  of  a  spatially  extended  scatterer  as  a 
fictitious  ensemble  of  identifiable,  stationary,  non-interactive ,  non- 
dispersive,  isotropic  point  seatterers).  It,  therefore,  avoids  the  problem 
of  the  conversion  of  radar  images  to  geometrical  images,  by  sidestepping  and 
avoiding  the  radar  image  altogether,  and  addressing  itself  to  the  problem  of 
generating  actual  three-dimensional  geometrical  images  directly  fron  radar 
data;  including  unique  optimal  results  from  incomplete  observation  aspect 
angles  and  frequency  information. 


A 

1 


I 


f 


* 


SECTION  I 


INTRODUCTION 


In  section  II,  the  general  inverse  scattering  and  radiation  problem 
associated  with  the  three-dimensional  inhomogeneous  scalar  field  Helmholtz 
wave  equation  is  formulated  as  a  Fredholm  integro-differential  equation  of 
the  second  kind  for  the  unknown  fields  and  sources  in  the  interior  of  a 
closed  surface  in  terms  of  the  known  'fields  on  this  closed  surface.  It  is 
shown  in  section  III  that  the  known  ansatz  to  this  integro-differential 
equation  reduces  to  an  angular  spectrum  Fourier  integral  if  the  closed  sur¬ 
face  integral  is  taken  over  the  far  fields  at  infinity.  In  section  IV,  this 
far  field  integro-differential  equation  is  reformulated  as  a  purely  algebraic 
equation  in  the  spatial  Fourier-transform  k-space,  yielding  a  purely  algebraic 
closed  form  solution  for  the  unknown  fields  in  this  k-space.  This  solution 
is  then  transformed  back  into  the  spatial  domain,  yielding  a  closed-form 
solution  to  the  far-field  integro-differential  equation,  consisting  of  a 
single  integral  resolvent  operator,  the  resolvent  kernel  of  which  is  the 
steady  state  Green’s  function  of  the  Poisson  equation  associated  with  the 
inhomogeneous  Helmholtz  equation.  This  scalar  inverse  theory  is  clearly 
applicable  to  quantum  mechanics,  acoustics,  etc. 

In  sect.  V,  the  inverse  integro-differential  equation  and  its  solu¬ 
tion  are  generalized  to  the  vector  wave  equation  resulting  from  Maxwell's 
equations . 

A  formal  physical  and  mathematical  definition  of  the  inverse  problem 
is  given  in  sect.  VI;  which,  in  essence,  defines  the  inverse  problem  as  that 
of  unknown  boundary  conditions  and/or  constitutive  equations  which  define 
the  specific  problem,  given  the  general  laws  and  the  specific  solutions 
which  are  the  observable  phenomenology. 
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The  uniqueness,  completeness,  and  well-behavedness  of  the  inverse 
solution  is  discussed  briefly  in  sect.  VII. 

It  is  shown  in  sect.  VIII  that  the  first  order  Neumann  series  solution 
of  the  inverse  integro-differential  equation,  as  well  as  the  first  order 
approximation  of  the  exact  solution  of  this  inverse  integro-differential 
equation,  is  identical  with  the  equations  governing  synthetic  microwave 
holography.  The  physical  interpretation  of  this  first  order  approximation 
is  discussed,  and  shown  to  be  equivalent  to  the  physical  optics  approximation. 

A  fast  Fourier  transform  method  for  evaluating  numerically  the  exact 
solution  to  the  inverse  problem  is  presented  in  sect.  IX. 

In  sect.  X,  the  three-dimensional  inverse  integro-differential  and  its 
closed  form  solution  are  generalized  to  four  dimensions,  with  the  surprisingly 
simple  result  of  the  resolvent  operator  being  algebraic.  This  four-dimension¬ 
al  solution  is  clearly  applicable  to  Doppler-contaminated  fields  scattered 
or  radiated  by  moving  scatterers. 

Since  a  distribution  theory  approach  was  taken  in  sect.  II  through  X, 
appendix  I  consists  of  a  simplified  and  unified  rederivation  of  the  classical 
Kirchhoff  solution  to  the  direct  radiation  problem,  presented  for  the  sake 
of  convincing  the  reader  of  the  validity,  elegance,  and  power  of  this  distri¬ 
bution  theory  approach. 

Appendix  II  consists  of  the  reformulation  of  Maxwell's  equations  into 
a  form  containing  only  the  fundamental  electric  and  magnetic  fields  (vis-a- 
vis  the  form  also  containing  the  displacement  and  induction  fields)  and  the 
total  charge  and  current  densities  (vis-a-vis  the  form  containing  the  free 
charge  and  current  densities  only),  thus  yielding  an  electromagnetic  wave 
equation  universally  valid  for  all  media.  Also  developed  in  this  appendix  II 
is  a  suitably  generalized  constitutive  equation  -  boundary  condition  relating 
this  total  current  density  to  those  fundamental  electromagnetic  fields. 

Both  these  reformulations  are  essential  for  the  subject  matter  covered  in 
sect.  V,  et.  seq. 
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SECTION  II 


I 


THE  INVERSE  SCATTERING  INTEGRAL  EQUATION 


Gauss'  theorem,  when  applied  to  the  fector  field  V(g4>),  yields 


dS*V(gsj>)  =  I  dv  V*V(g<J>) 


dv  V2(g<ji)  , 


provided  VCg<J>)  is  continuous  and,  hence,  differentiable  in  v  and  on  s.  These 
continuity  conditions  can  be  totally  dispensed  with  if  the  scalar  fields  $ 
and  g,  and,  hence,  V(g<{>),  are  taken  as  distributions  [17];  specifically,  if 
<j>  is  taken  as  a  field  satisfying  the  inhomogeneous  (neither  source-  nor 
singularity-free)  Helmholtz  (time-reduced  wave)  equation 


V2(j>  +  ko  $  =  -  p  , 


and  g.  is  taken  as  the  associated  free-space  Green's  function 


S<x|x'>  =  , 


r  =  x-x'  , 


r  =  r 


mr- 


which  satisfies  the  inhomogeneous  differential  equation 


V2g  +  ki  g  =  -  6  f 


(7) 


and  where  both  and  g  satisfy  the  radiation  condition  at  infinity,  (The 
free-space  wave  number  — -  =  ^  is  designated  by  k©,  vis-a-vis  the  conventional 
notation  k,  since  the  latter  will  subsequently  be  used  to  designate  the 
Fourier  transform  variable  of  the  spatial  coordinate  x) . 

In  order  to  further  convince  the  reader  not  familiar  with  distribution 
theory  and  its  validity  and  power  when  applied  to  field  problems,  this  same 
distribution  theory  approach  is  applied  in  appendix  I  to. the  rederivation  of 
the  classical  Kirchhoff  method  of  integration  of  the  field  equations  for 
the  direct  radiation  (and  scattering)  problem,  without  imposing  the  continuity 
restrictions,  and,  hence,  without  the  classical  singularity  isolation  spheri¬ 
cal  surface  which  is  taken  to  the  limit  of  vanishing  size.  The  derivation  of 
the  inverse  scattering  integral  equation  presented  in  this  sect.  II  can  also 
be  accomplished  by  such  classical  means;  however,  the  resulting  mathematical 
derivation  is  beset  with  vastly  increased  cumbersome  details,  which  only 
obscure  the  physical  meanings  involved.  ;■ 

The  left-hand  side  of  (1)  yields 


^ dS*V(g$)  =  ^ 


dS*VCg$)  =  <f>  (g  +  4  Vg)*dS 
s  s 


=  f(s!£+  ’ 


(8) 

(9) 


and  the  right-hand  side  of  (1)  yields 


dv  V* V<g«J>)  =  dv  V*(g  V4>  +  <{>  Vg) 
=  j  dv  [V 


(g  7<(>)  +  V •  C 4>  Vg)] 


Cio) 


(ii) 


=  I  dv  <Vg-V<{>  +  g  V2<j>  +  V<|>*Vg  +  $  V2g)  (12) 

*v 

=  j  dv  (g  V24>  +  <j>  V2g  +  2V$‘Vg)  ,  (13) 

•'v 


which,  with  the  aid  of  (3)  and  (7)  yields 


J  dv  V2(g$)  =  J  dv  [g(-kl$-p)  +  <j>(-k!g-S)  +  2Vg-V<j>] 

=  I  dv  (-klg<|>-gp-kog<j>-6<i>+2Vg*v<j>) 

^  v 

=  j  dv  C  2Vg  *  V4>  -  2k|  g<{>-gp-S<j>) 

J  v 

=  I  dv(2Vg*Vi}>-2k£g<j>)  -  j  dv  gp  -  j  dv  6<|> 

*  v  V  v 


(14) 


(15) 


(16) 


(17) 


If  the  field-point  X*  is  taken  as  inside  the  volume  bound  by  the 
surface  s,  then 


I 


dv  6  4>  = 


(18) 


By  the  conventional  Kirchhoff  interpretation  of  (3)  and  (7)  for  the 
internal  (to  s)  field  point  X'  (see  appendix  I), 
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(19) 


I 


f 


dv  g  p  =  4>*  +  Vg  —  g  V$)*dS 


where  the  right-hand  side  Kirchhoff  surface  integral  of  (13)  represents  the 
contribution  to  the  field  at  the  field  point  X'  inside  s  due  to  sources  p 
outside  of  s. 

Thus,  if  all  the  sources  are  taken  as  in  v  (i.e.,  inside  s;  i.e., 
p (X)=0  for  all  X  outside  s) ,  then  (19)  yields 


dv  g  p  =  <j>’ 


(20) 


By  a  similar  argument,  if  the  sources  outside  of  s  do  not  vanish,  then 
these  sources  external  to  s  must  be  taken  as  giving  rise  to  an  externally 
imposed  (external  to  the  inverse  problem)  incident  field  <{>’.  ,  given  by  (19) 

S' 

as 


^  (<f>  Vg  -  g  )  »dS 
s 


(21) 


Equation  (19)  thus  yields 

(20.1) 

(20.2) 

i 

where  6  is  the  scattered  field  at  the  field  point  X*.  Since  the  scattered 

c 

field  also  satisfies  the  Helmholtz  equation  (3),  equations  (1)  through  (20) 
remain  valid  if  the  field  is  taken  as  the  scattered  field  *  only.  Thus 
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no  formal  distinction  arises  between  the  inverse  radiation  (active)  and 
inverse  scattering  (passive)  problems. 


With  the  aid  of  (18)  and  (20),  (17)  reduces  to 


dv  (Vg*V$>  -  k§  g  <f>)  -  2<{> * 

v  v 

With  the  aid  of  (8)  and  (22) ,  (2)  reduces  to 


(22) 


dv  (Vg*V<|)  -  k|  g 


-f 


(g  Vij>  +  (ji  vg) *ds 


(23) 


In  the  context  of  the  inverse  problem,  the  integrand  of  the  right-hand 
side  surface  integral  of  (23)  is  measurable,  and,  hence,  this  surface  integral 
is  evaluable;  it  thus  becomes  convenient  to  take  this  right-hand  side  of  (23) 
as  the  known  ansatz  e  to  the  inverse  problem;  i.e.. 


0>  =  -  1 


By  (23)  and  (24),  the  inverse  problem  can  thus  be  stated  by  the 
Fredholm  integro-differential  equation  of  the  second  kind 


f 


(g  V<J>  +  <}>  Vg)»dS 


(24) 


(25) 


It  should  be  noted  that  the  classical  Kirchhoff  surface  integral  (see 
appendix  I)  vanishes  in  general  if  both  the  field  point  as  well  as  all  the 
sources  are  inside  the  surface  (i.e.,  X*€  s  and  p(X)=0  for  all  X  ^  s),  where 
as  the  ansatz  surface  integral  (24)  does  not  vanish  in  general  for  these  same 


conditions.  The  physical  reason  for  the  vanishing  of  the  Kirchhoff  surface 
integral  is  clearly  that  there  are  no  sources  outside  of  this  surface;  since 
this  surface  integral  represents  the  fields  produced  by  sources  outside  of 
this  surface.  The  mathematical  reason  for  this  vanishing  is  simply  that  the 
(-)  sign  in  the  integrand  renders  the  two  integral  terms  as  canceling  each 
other.  Whereas  the  ansatz  surface  integral  does  not  Vanish  since  these  same 
two  terms  add  up  by  virtue  of  the  (+)  sign  in  the  integrand  (if  the  ansatz 
is  in  the  far  field,  then  these  two  integrands  (vis-a-vis  integrals)  are 
equal).  It  is  indeed  remarkable  that  the  sole  distinction  between  these 
two  surface  integrals  should  turn  out  to  be  the  (+)  vs.  (-)  sign  in  the 
integrand,  and  that  this  sole  transition  from  (-)  to  (+)  should  render  the 
Kirchhoff  surface  integral  as  a  useful  ansatz  (surface  integral)  to  the 
inverse  problem. 


SECTION  III 


THE  FAR-FIELD  INVERSE  SCATTERING  INTEGRAL  EQUATION 


If  only  the  far- field  is  known,  then  the  surface  integral  of  the 
ansatz  (24)  can  clearly  be  taken  as  a  spherical  surface  at  infinity.  Let 
R  be  the  radius  of  this  surface  (see  fig.  11);  it  thus  follows  that 

e’  =  -  {  r  (g  V4>  +  d>  vg) *R  r2  dn  ,  (26) 

i 

where  the  implied  integration  is  now  over  the  complete  solid  angle  0. 

The  gradient  of  the  Green's  function  (see  eqn.  (4))  is 

ik.r 

Vg  =  (ik8r  -  D  r  ,  (27) 


which,  in  the  far-field  (retaining  only  terms  in  ~)  reduces  to 


-  J-k0r 
e 


vg  =  ik.r  ^ 


=  iko 


3ik0r 

4nr 


(28) 


(29) 


where  k#  is  the  far-field  observation  propagation  wave-number  vector  (see 
fig.  11). 


FIGURE  11 

FAR- FI  ELD  GEOMETRY 


The  field  4>  in  the  far-field  (retaining  only  terms  in  — ) ,  can  be 

R 

expressed  by  virtue  of  the  radiation  condition  at  infinity,  as  (see  fig.  11) 


$(x)  -  <l>(k0)  , 

R 


where  ^Ck»)  is  now  the  range  and  phase  normalized  (relative  to  a  given 
coordinate  system  in  which  the  far-field  is  measured  and  the  near  field  <j>  is 
to  be  determined;  see  fig.  11).  \J>(k0)  clearly  depends  on  the  observation 

direction  ko  only.  The  factor  ^4tT  was  chosen  for  the  sake  of  consistency  with 
the  conventional  definition  of  power  cross  section. 


The  gradient  of  this  far-field  (30)  is  thus 


_ikoR  ik0R 

e  ..  a 

=  -  (ikoP  -1)R  t}i(ko )  + - 


7r|»(k0) 


which,  in  the  f&r-field  (retaining  only  terms  in  reduces  to 

K 


eikoR 

v*  - -  [ikoR  ij>(ke)  +  V'Mko)] 

/4tT  R 
eikoR 

[iko^(ko>  +  vij»(ko)3 


fa  R 


(32) 


(33) 


With  the  aid  of  (4) ,  (29) ,  (30) ,  and  (33) ,  (26)  reduces  to 


/*  ik  r  e*koR 

6f  -  -  *  R^!j  [2iko*(ko)+v*(ko>]«R  R2  do 


(34) 


=  -  i  (47rr3/2  L*fl  f  eiko(r+R)  B  [2ik»*R  *(ko>+R-'V*(k6-)3  dO  .  (35) 


Jr, 


Since  ^(k0)  varies  only  in  the  angular  (observation)  direction,  and 
not  in  the  radial  direction,  it  follows  that  (see  fig.  11) 


R*Vtj;(ko>  =  0 


(36) 


Also  (see  fig.  11) 


k0-R 


ko 


(37) 


Equation  (35)  thus  reduces  to 


6’  =  -  iko  (4tt)~3/2  ^mit  f 9iko(r+R)  B  ^,(ko) 

K  +  j  r 

•'ll 


dO 


(38) 
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Examination  of  fig.  11  readily  reveals  that 


Limit  R  ik0(r+R)  2ik0R  -ik0*X* 
—  e  —  e  e 


R  +  ®  r 


Equation  (38)  thus  reduces  to 


8’  =  -  ik0  (4tt)“3/2  e2lk 


,r  j  Q-iko*x' 


^(k0 >  dn 


Since  the  arbitrary  phase  factor  of  (2k0R)  is  common  to  all  the 
measured  ansatz,  it  becomes  convenient  to  arbitrarily  chose  it  as  an  integer 
multiple  of  2ir,  thereby  reducing  (40)  to 


0’  *  -  ik0  (4ir) 


~3/ 2  f  -iko-X' 


Tp(ko)  dn 


It  further  becomes  convenient  to  now  interchange  the  variables  x  and 
X'  in  (25)  and  (41);  thus  yielding  the  far-field  integro-differential  equation 
and  its  ansatz  respectively  as 


where  the  volume  integral  in  (42),  bound  by  the  surface  at  infinity,  is  now 
over  all  space. 
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SECTION  IV 


SOLUTION  0?  THE  FAR-FIELD  INVERSE  SCATTERING  INTEGRAL  EQUATION 


In  cartesian  coordinates,  the  Green’s  function  g(x|x')  can  be  written 
as  (see  (4),  (5),  and  (6)) 


g(x}x? )  =  g(X-X’) 


(44) 


The  integro-differential  equation  (42)  can  thus  be  written  as  the 
three-dimensional  convolution  equation 


I 


<j>(X)  -  III  [V’gCX-X’^V'^X’)  -  k!  g(X-X')  <j>(X’)]  d3*'  -  0(X)  ,  (45) 


which,  in  a  more  compact  notation,  is 


4  -  Vg?V<j>  +  k!  g#ij>  =  6 


(46) 


where  the  operation  (*)  designates  the  vector  dot-product  convolution. 

If  the  fields  and  the  Green's  function  in  (46)  are  taken  as  distri¬ 
butions  for  the  sake  of  the  existence  of  their  three-dimensional  Fourier 
transform;  i.e.. 
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' v 

’  ’  V 

Vv  ' 

■  1 

H.  i 

•  * 

■ 

‘  ‘V  • 

1 

V  ‘  ■ 

> 

B  *(X)  +■*■  $>(k)  (47)  j 

■  8(X)  ++  0(k)  (48) 

H  g(X)  -M-  G(k)  ,  (49) 

%  *  *v 

H|  and,  by  the  three-dimensional  differentiation  rule  for  Fourier  transforms 

■  V+(X)  <->  ik  $(k)  (50) 

'•  ['  ■ 

• . 

B  Vg(X)  -w  ik  G(k)  ,  (51) 

'  ,  *  j  > 

B 

'  J 

Bl  then,  by  the  three-dimensional  convolution  theorem,  the  three-dimensional 

c 

m\  Fourier  transform  of  (46)  yields  in  the  three-dimensional  transform  k-space 

'A'";  , •* 

raj  the  algebraic  equation 

V  • 

♦ 

,  .  -# 

B  *(k)  +  k2  G(k)  *(k)  +  kl  G(k)  *(k)  =  o(k)  (52) 

■  [1  +  (k2  +  k§)  G  ]  4  =  0  .  (53) 

«  .  ** 

■  .# 

raj >  Furthermore,  since  for  the  three-dimensional  delta-function  taken  as  ' 

Bl  a  distribution  •; 

■  6U5  -w  1  ,  (54) 

B||  it  follows  that  the  three-dimensional  Fourier  transform  of  (7)  yields  in 

H|  k-space  the  algebraic  equation  ^  - 

•*  .  '  'S’ 

'  •  *  # ' 

-----  -  “r  ■■ . '  — | 
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-  k2  G(k>  +  kl  G(k)  =  -  1 


G(k)  = 


k2  -  kl 


With  the  aid  of  (56),  (53)  yields 


$(k)  +  (k2  +  kl)  GCk)  4>(k)  =  0(k) 


(k2  +  kl  \ 

- )  *(k>  =  0(k) 

k2  -  kl/ 


k2  +  kl\ 

1  + - j  =  G(k) 

k2  -  kl/ 

k2  -  kl  k2  +  kl  \ 

- + - )  $(k)  =  o(k) 

k2  -  kl  k2  -  kl  / 


k2  -  kl 


$<k)  =  e(k) 


/  k2  -  kl  \ 

$(k)  =  ( - )  ©ck) 

\  2  k2  / 


=  £  ©(k)  -  ~  0(k)  . 


Next,  the  steady-state  Green's  function  g«,  given  by 


g«(x) 


1 

4irr 


(64) 


satisfies  the  Poisson  differential  equation  (associated  with  the  Helmholtz 
equation  (3)) 


V2g0(X)  =  -  S(X) 


(65) 


Again,  taking  this  steady-state  Green's  function  as  a  distribution, 
assures  the  existence  of  the  three-dimensional  Fourier  transform  of  this 
steady-state  Green's  function;  i.e.. 


g0(x)  G„(k) 


(66) 


It  thus  follows  from  (65)  by  the  three-dimensional  differentiation 
rule  that 


-  k2G6(k)  =  -  1 


(67) 


Go(k)  =  £2 


(68) 


Equation  (63)  can  thus  be  written  as 


k§ 


»(k)  =  i  ©(k)  -  f 


Go (k)  ©(k) 


(69) 


58 


and  by  the  three-dimensional  convolution  theorem,  it  thus  follows  that  (still 
in  cartesian  coordinates  only) 


<k(X)  =  i  e(X)  -  g0(X)*6(X) 


(70) 


I 


=  i  0(x)  -  ^  III  go(X-X*)  e(x’)  d  x* 


(71) 


and,  consequently,  in  any  coordinate  system 


$(X)  =  *  6(X) 


-  f  f  g.<x|x* 


)  e(x’)  dvf  . 


(72) 


It  thus  follows  by  (64)  that  the  solution  to  the  integro-differential 
equation  (42)  is 


(73) 


It  is  remarkable  that  the  resolvent  kernel  [18]  of  the  integro- 
differential  equation  (42)  should  turn  out  to  be  the  steady-state  Green's 
function;  the  physical  implications  of  this  fact,  which  are  yet  unclear, 
should  be  explored. 


SECTION  V 


THE  GENERALIZATION  FOR  THE  ELECTROMAGNETIC  VECTOR  FIELDS 


The  generalization  of  the  previously  developed  scalar  field  inverse 
scattering  theory  to  a  vector  field  inverse  scattering  theory  applicable  to 
the  electromagnetic  field  can  be  accomplished  analogously  to  the  generaliza¬ 
tion  of  Stratton  and  Chu  [19]  and  Franz  [20]  of  the  Kirchhoff  method  of  inte¬ 
gration  of  the  scalar  field  equations  to  vector  fields,  applicable  to  direct 
scattering  or  radiation;  i.e.,  by  replacing  (1)  with  Gauss'  theorem  of  the 
form 

d$xV*(gE) 

dSxVx(gH) 

where  the  electric  and  magnetic  vector  fields  E  and  H  satisfying  the  vector 
wave  equations 


■I 


dv  VxVx(gE) 


(74) 


•/ 


dv  VxVx(gH) 


(75) 


VxVxE  -  ko  E  =  iupo  d  (76) 

VxVxH  -  kl  H  =  vxd  ,  (77) 


which  replace  the  scalar  wave  equation  (3)  used  previously. 
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From  here  on,  the  generalization  proceeds  in  a  fashion  totally  analogous 
to  the  previously  presented  development;  however,  with  the  following  one 
important  exception.  The  scalar  theory  developed  consists  of  a  single  field, 
and,  hence,  a  single  statement  of  Gauss'  theorem  (1)  and  a  single  wave  equation 
09 .  The  vector  theory  consists  of  two  fields  (the  electric  and  magnetic 
fields),  and,  hence,  a  dual  statement  of  Gauss'  theorem  (74),  (75)  and  two 
wave  equations  (76),  (77).  This  dual  set  of  vector  equations,  is,  however, 
coupled  by  a  common  source  (the  current  density  J  in  the  wave  equations  (76) , 
(77)),  as  well  as  Maxwell's  equations,  which  are  a  sot  of  first-order 
differential  equations  for  which  there  is  no  scalar  analog,  and  which  are 
used  to  entwine  the  coupling. 

Regretably,  however,  such  a  procedure,  although  executable,  tends  to 
bog  down  in  massive,  cumbersome,  and  inelegant  details,  which  tend  to  obscure 
the  physical  meanings  involved.  An  alternate,  much  simpler  derivation,  which 
yields  the  identical  results,  will  thus  be  presented  next. 

Prior  to  proceeding  with  this  derivation,  however,  a  digressionary 
discussion  of  the  form  of  Maxwell's  equation  most  appropriate  to  the  inverse 
problem  is  in  order.  For  organizational  reasons,  however,  this  discussion 
is  relegated  to  appendix  II,  to  which  the  reader  is  referred. 

The  universally  valid  time-reduced  vector  wave  equation  for  the 
magnetic  vector  field  H,  in  terms  of  the  total  currents  J  (see  appendix  II)  is 


VxVxH  -  k§  H  *  VxJ 


(78) 


which  reduces  to 


w-H  -  v2H  -  kl  H  =  vxj  , 


(79) 
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which,  by  Maxwell's  second  equation  (see  appendix  II,  eqn.  (33)),  further 
reduces  to 


V2H  +  k§  H  =  -  VxJ 


(80) 


If  a  source  vector  p  is  defined  as 


P  5  VxJ  , 


(81) 


then  (80)  can  be  written  as 


V2H  +  k§  H  =  -  p 


(82-) 


If  and  designate  the  cartesian  components  of  the  vectors  H  and 
P  respectively  (where  y=1 ,2,3.) ,  then  (82)  can  be  rewritten  in  a  mixed  vector  - 
cartesian  tensor  notation  as 


V2Hp  +  kl  Hv  =  -  Pjj  ;  y=1,2,3.  (83) 

Each  of  the  three  equations  (83)  clearly  is  of  the  form  of  the  scalar 
field  wave  equation  (3) ,  it  thus  follows  that  the  entire  previously  developed 
scalar  field  inverse  scattering  theory  is  applicable  to  each  of  the  three 
cartesian  components  of  H  and  p .  The  inverse  scattering  integro-differential 
equation  (25)  and  its  ansatz  (24)  can  thus  be  written  (still  in  mixed  vector  - 
cartesian  tensor  notation)  respectively  as 
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I 


H'  -  I  dv  (Vg*7H  -  k|  g  H  )  =  e'  ,  y= 1 ,2,3. 

v  J  yP  ay  v 


dS-(g  VHy  +  Hp  Vg)  ,  y=1,2,3. 


which  clearly  can  be  reconciled  into  the  pure  vector  notation  form 


H’  -  J  dv  (Vg-VH  -  k|  g  H)  =  6' 


-f 


dS*  (g  vH  +  Vg  H) 


Similarly,  the  far-field  inverse  scattering  integral  equation  (42)  and 
its  ansatz  (43)  yield  respectively 


H  -  I  dv’  (v'g'.v'H*  -  k§  g'  H’)  =  0 


-  ik0  (4tt)~3/2  j  e“ik°*X  ^(ko)  df 


where,  by  (30),  the  range  and  phase  normalized  magnetic  far-field  $(ko)  is 
given  by 


H(X)  =  -  \Kk0> 

v^T  R 


and  again  similarly,  and  finally,  the  solution  (73)  to  this  far-field  inverse 
scattering  integral  equation  yields 
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(91) 


Unfortunately,  however,  no  such  simple  derivation  for  the  electric 
field  has  been  found  to  date.  However,  commencing  with  (74),  (76),  et  seq., 
a  set  of  electric  field  equations  totally  analogous  to  (86)  through  (91)  is 
derivable;  i.e.,  a  set  of  equations  (86)  through  (91)  in  which  the  magnetic 
field  H  is  replaced  with  the  electTic  field  E.  The  details  of  this 
derivation  are  not  presented  here  because  of  their  cumbersomeness,  and 
because  these  details  shed  no  light  on  the  physical  meanings  involved. 
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SECTION  VI 


FORMAL  STATEMENT  AND  SOLUTION  OF  THE  INVERSE  PROBLEM 


An  academic  digression  to  a  formal  statement  and  definition  of  the 
inverse  problem  is  undertaken  next.  First,  however,  a  careful  examination, 
formal  restatement,  and  definition  of  the  classical  direct  problem  is  in 
order. 


From  a  mathematical  point  of  view,  a  direct  problem  consists  of  a 
given  differential  equation,  a  given  set  of  boundary  conditions,  and  an 
unknown  solution  (or  spectrum  of  solutions).  From  a  physical  point  of  view, 
the  differential  equation  represents  the  universally  valid  law,  the  boundary 
conditions  represent  the  nature  of  the  specific  problem,  and  the  unknowns  are 
the  observable  phenomenology.  Additional  insight  is  gained  by  examining  the 
structure  of  the  equivalent  integral  equation  formulation  of  the  direct 
problem.  Mathematically,  such  integral  equation  formulations  are  in  general 
arrived  at  in  the  following  fashion.  An  integral  representation  of  the 
differential  equation  is  obtained  (usually  by  an  integration  of  this 
differential  equation  with  the  aid  of  an  appropriate  Green's  function),  which 
is  then  combined  with  the  appropriate  boundary  conditions  to  yield  an  integral 
equation.  Physically,  the  integral  representation  of  the  differential  equation 
is  merely  an  integral  restatement  of  the  universally  valid  law,  which,  in 
general,  clearly  illustrates  its  global  nature  (i.e.,  that  fields  at  one 
point  are  determined  by  sources  at  all  other  points).  Conversely,  however, 
the  boundary  conditions  are  a  restatement  of  the  specific  constitutive 
equations,  which  determine  the  specific  problem.  Such  constitutive  equations 
are  in  general  local  in  nature  (i.e.,  the  specific  relation  between  the 
sources  at  one  point  and  the  fields  at  that  same  point  only) .  The  following 
equivalence  between  boundary  conditions  and  constitutive  equations  is  thus 
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of  importance.  Mathematically,  a  constitutive  equation  is  a  spatially 
distributed  boundary  condition;  and,  conversely,  physically,  a  boundary 
condition  is  a  localized  constitutive  equation.  (Such  boundary  conditions  - 
constitutive  equations  in  general  arise  from  a  localized  imposition  of  the 
appropriate  laws  of  motion  governing  the  source  distribution.)  From  an 
integral  equation  point  of  view,  it  is  this  boundary  condition  -  constitutive 
equation  that  determines  and  imposes  the  given  known  domain  of  integration, 
and,  hence,  the  nature  of  the  specific  problem. 

From  a  mathematical  point  of  view,  an  inverse  problem  thus  consists  of 
a  given  differential  equation  or  its  equivalent  integral  representation,  a 
given  solution  (or  spectrum  of  solutions) ,  and  an  unknown  boundary  condition 
(or  its  equivalent  domain  of  integration) .  From  a  physical  point  of  view, 
an  inverse  problem  thus  consists  of  a  given  universal  law,  a  given  (observed, 
measured)  phenomenology,  and  an  unknown  constitutive  equation. 

Having  defined  the  direct  and  inverse  problems  in  the  preceding 
fashion,  it  follows  naturally  that  if  the  differential  equation  (or  its 
integral  representation)  stating  the  universal  law  is  the  unknown,  then  the 
problem  can  be  classified  as  the  basio  or  fundamental  problem  (there  exists 
no  formal  methodology  to  date  for  solving  such  problems) . 

The  following  short  table  is  thus  in  order: 


PROBLEM  CLASSIFICATION 


MATHEMATICAL  ENTITY 

PHYSICAL  MEANING 

IF  UNKNOWN,  PROBLEM  IS 

Differential  Equation 
(or  equivalent  integral 
equation  representation) 

Universal  (Global)  Law 

Basic,  Fundamental 

Boundary  Condition  (or 
domain  of  integration) 

Constitutive  (Local) 
Equation 

Inverse 

Specific  Solution  to 
differential  (or 
integral)  equation 

Observable  (Measurable) 
Phenomenology 

Direct 
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The  inverse  scattering  (or  radiation)  problem  is  thus  the  problem  of 
determining  the  unknown  boundary  conditions  -  constitutive  equations  (i.e., 
the  structure  of  the  scattered  or  radiator) ,  given  the  field  equations  and 
the  incident  and  scattered  (or  tot?*  radiated)  fields.  The  problem  of  deter¬ 
mining  the  unknown  source  distribution  that  gave  rise  to  the  observed 
scattered  field  (or  radiation),  although  intimately  connected  with  the 
inverse  scattering  (or  radiation)  problem,  is  thus  only  incidental  to,  and 
not  the  final  objective  of,  the  inverse  scattering  (or  radiation)  problem. 

The  distinction  between  the  inverse  scattering  and  inverse  radiation  problem 
is  the  mathematical  and  physical  triviality  of  whether  the  unknown  structure 
is  radiating  passively  or  actively,  respectively  (i.e.,  whether  the  observed 
fields  are  measured  actively  or  passively,  respectively) . 

For  the  direct  vs.  the  inverse  scattering  (and  radiation)  problem,  the 
following  specific  additional  mathematically  fundamental  topological 
distinction  arises.  For  the  direct  problem,  the  Kirchhoff  closed  surface 
integral  ^>ds*(g  V<j>  —  4>  Vg)  relates  fields  at  a  field  point  located  on  one 
side  of  the  closed  surface  to  all  the  source  distribution  located  on  the  other 
side  of  this  closed  surface  (see  appendix  I).  For  precisely  this  reason, 
since  the  inverse  problem  is  characterized  by  unknown  fields  as  well  as 
unknown  sources,  this  Kirchhoff  surface  integral  is  thus  rendered  totally 
useless  to  the  inverse  problem.  However,  the  closed  surface  integral 
^dS*(g  V4>  +  7g)  relates  the  fields  at  a  field  point  located  on  one  side  of 

the  closed  surface  to  the  source  distribution  located  on  the  same  side  of 
this  closed  surface  (see  (24)  and  (25)).  It  is  for  this  reason  that  this 
latter  closed  surface  integral  is  rendered  a  useful  ansatz  to  the  inverse 
problem. 

In  the  context  of  this  definition  of  the  inverse  problem,  it  is  thus 
clear  that  neither  a  solution  for  the  electric  and  magnetic  fields  (see  (86) 
through  (91)),  nor  a  solution  for  the  total  current  distribution  J  that  gave 
rise  to  these  fields,  is  a  solution  to  the  inverse  problem;  however,  a 
solution  for  the  unknown  boundary  conditions  -  constitutive  equations  which 
locally  connect  these  electric  and  magnetic  fields  to  the  current  distri¬ 
bution,  is  a  solution  to  the  inverse  problem.  It  is  to  this  latter  problem 
that  the  remainder  of  this  section  is  addressed. 
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Given  a  solution  for  E  and  H  by  (86)  through  (91) ,  the  total  current 
distribution  J  can  be  computed  on  a  local  basis  via  the  universally  valid 
Maxwell's  fourth  equation  (see  appendix  II,  eqn.  (35)  and  (46))  for  the  total 
current ;  i . e . , 

d  =  -  itdeo  E  +  VxH  .  (92) 

Thus,  given  the  electric  and  magnetic  fields,  as  well  as  the  total 
current  density,  the  coefficients  of  the  generalized  constitutive  equation 
(see  appendix  II,  eqn,  (49)),  which  determine  (reveal)  the  structure, 
geometry,  and  electromagnetic  properties  of  the  scatterer  (or  radiator) ,  can 
be  readily  computed.  In  fact,  since  these  scalar  constitutive  coefficients 
in  the  vector  constitutive  equation  are  over-determined  by  knowledge  of  the 
vector  fields,  a  least-square-best-estimate  (in  the  presence  of  noise)  of 
these  coefficients  can  be  determined  as  well. 


SECTION  VII 


UNIQUENESS,  COMPLETENESS,  AND  WELLBEHAVEDNESS  OF  THE  INVERSE  SOLUTION 


The  problem  of  uniqueness,  completeness,  and  wellbehavedness  of  the 
inverse  solution  (under  incomplete  and  noisy  error  contaminated  ansatz)  is 
intimately  connected  with  the  algebraic  nature  of  the  k-space  transform 
representations  (57)  and  (69)  of  the  inverse  integral  equation  and  its 
solution  reprectively,  the  Fourier  surface  integral  transform  nature  of  the 
associated  ansatz  (43) ,  and  the  extension  of  the  entire  inverse  theory  and 
its  solution  to  the  wide-band  (short -pulse)  case.  For  these  reasons,  both 
the  method  of  derivation  as  well  as  the  final  results  governing  the  uniqueness, 
completeness,  and  wellbehavedness  of  the  solution  are  identical  to  the  method 
of  derivation  and  results  governing  the  physical  optics  inverse  scattering 
theory  (see  part  I),  A  further  detailed  investigation  and  study  of  this 
subject  is,  however,  currently  being  undertaken  under  a  follow-on  contract 
(to  contract  covered  by  this  report).  A  detailed  treatment  of  this  subject 
is  thus  being  deferred  to  subsequent  report (s)  covering  this  follow-on 
contract . 
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SECTION  VIII 


> 


THE  FIRST-ORDER  APPROXIMATION. 


The  first-order  approximation  Neumann  series  solution  of  the  general 
inverse  integro-differential  equation  (23) ,  with  the  aid  of  (24) ,  clearly  is 


+  '  = 


(g  V£  +  Vg)*d$ 


(93) 


which,  for  the  far-field  ansatz  case,  with  the  aid  of  (43),  reduces  to 


*(X)  =  -  ik0  (4tt)"3/2  j  e'iko*X  t|>(ke)  d ft  ;  (94) 

which  turns  out  to  be  identical  to  the  first-order  approximation  of  the  exact 
far-field  solution  (73) . 

This  first-order  solution  is  identical  to  the  synthetic  microitaoe 
holography  solution  of  Tricoles  [21],  which  has  already  enjoyed  considerable 
experimental  verification.  The  derivation  of  this  first-order  solution  can 
thus  be  taken  as  the  formal  theoretical  foundation  of  this  synthetic  micro- 
wave  holography. 

In  the  derivation  of  the  synthetic  microwave  holography  equations,  the 
’’ansatz*'  is  taken  as  (£ds*(g  V<?*  -  Vg)  where  <{>*  denotes  the  complex 
conjugate  of  the  outwardly  radiated  field;  which,  in  turn,  heuristically 
represents  the  "direction  of  propagation  reversed"  equivalent  "inwardly 
focussed"  field.  From  a  holographic  point  of  view,  this  ansatz  also  could 
have  been  taken  as  (£  dS*Cg*  V4>  —  Vg*),  where  g*  denotes  the  complex 
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conjugate  of  the  ••outgoing”  Green's  function;  which  in  turn  represents  the 
"ingoing”  Green's  function.  Since  the  complex  conjugate  of  both  the  field 
as  well  as  the  Green’s  function  satisfy  differential  equations  similar  to  (3) 
and  (7)  respectively,  i.e.. 


+  k!  <f>*  =  -  p*  (95) 

V2g*  +  k§  g*  =  -  6  ,  (96) 


it  thus  follows  that  the  entire  development  of  this  report  could  also  have 
been  done  in  terms  of  either  this  complex  conjugate  field  or  this  complex 
conjugate  Green's  function;  with  the  sole  exception  that  the  ansatz  would 
have  turned  out  to  be  the  above  holographic  ansatz  with  the  (-)  sign  in  the 
integrand  (vis-a-vis  the  (+)  sign  in  the  integrand  of  (93)).  It  should  be 
pointed  out  that  in  the  case  of  the  electromagnetic  vector  field  ansatz  (74) 
and  (75),  only  conjugation  of  the  Green's  function  yields  the  current  results. 
Such  an  approach,  although  heuristically  somewhat  attractive,  tends,  however, 
to  burden  the  rigorous  derivation  of  the  exact  inverse  formulation  and  its 
solution  with  mathematically  cumbersome  details,  which  also  obscure  the 
physical  meanings  involved. 

An  investigation  of  the  properties,  limitations,  and  physical  meaning 
of  this  first-order  approximation  (93)  and  (94)  will  be  taken  up  next.  To 
this  end,  let  the  source  distribution  p  be  computed  from  this  first-order 
approximation;  i.e.,  by  the  basic  wave  equation  (3), 


p  ~  *  v2$  -  k!  $ 


(97) 


thus,  by  (94) 
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P  =  iko  (4tt)"3/2  (V2  +  kl)  |  e  ij»(k0)  dn 

n 


=  iko  (4tt)"3/2  f  (v2  +  ki)  e"iko*x  g>(ke)  dn 

n 

=  iko  (4tt)"3/2  J* (-kl  +  kl)  e“iko*X  4>(k0)  dQ 


=  0 


(98) 


(99) 


(100) 


(101) 


The  same  result  of  (101)  can  be  obtained,  although  mathematically  much 
more  laboriously,  by  replacing  (94)  with  (93)  in  (98). 

It  thus  follows  from  (101)  that  the  source  distribution  is  not 
reconstructable  from  the  first-order  solution;  and  that,  the  fields  reconstructed 
by  this  first-order  solution  satisfy  the  homogeneous  free-space  wave  equation 


V2$  +  kl  <J>  =  0 


(102) 


This  first-order  reconstructed  source-free  field,  satisfying  the  free- 
space  wave  equation,  is  thus  consistent  with  synthetic  holography  interpreta¬ 
tion  of  the  first-order  solution  (93) .  The  physical  meaning  of  this  first- 
order  solution  is  thus  that  of  first-order  scattering  (or  radiation) ,  and 
consistent  with  the  physical  optics  approximation  (direct  !,rays"  only) . 

The  reconstruction  of  the  unknown  boundary  condition  -  constitutive 
equation  (see  sect.  VI)  which  determines  the  unknown  geometrical  structure 
and  unknown  electromagnetic  properties  of  the  scatterer  (or  radiator)  requires 
knowledge  of  both  the  fields  and  the  sources.  Thus,  since  the  first-order 
solution  cannot  reconstruct  the  sources,  it  follows  that  from  an  inverse 
scattering  point  of  view  (see  sect,  VI),  that  the  fiTSt-order  solution 
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possesses  the  critical  pathology  of  not  being  an  inverse  solution  at  all. 

Examination  of  the  complete  solution  (73)  thus  reveals  that  all  the 
source  distribution  information,  and,  hence,  all  the  formal  inverse  scattering 
information,  is  contained  only  in  the  resolvent  operator  of  (73);  i.e.,  in 
all  the  higher- than- first-order  Neumann  series  solution  terms. 


SECTION  IX 


NUMERICAL  EVALUATION  OF  THE  RESOLVENT  OPERATOR 
BY  MEANS  OF  THE  FAST  FOURIER  TRANSFORM 


Let  N  be  the  number  of  data  points  in  the  three-dimensional  space  x 
for  which  the  ansatz  6(X)  has  been  evaluated  by  (43).  Evaluation  of  the  field 
<}>(X)  for  these  N  points  in  X  by  conventional  numerical  integration,  as  implied 
by  the  resolvent  operator  in  (73),  will  thus  require  N2  multiply-add  operations 
and  storage  allocations.  To  recognize  that  this  is  an  untenable  situation 
requires  merely  an  examination  of  the  practical  order  of  magnitude  of  N. 

A  minimally  reasonable  three-dimensional  resolution  of  the  field  is  of  the 
order  of  100  points  in  each  of  the  three  dimensions';  i.e.,  N= 1 06  and  N2=1012, 
which  is  totally  unacceptable  in  practice  for  even  the  largest  and  fastest 
state-of-the-art  computers. 

In  cartesian  coordinates,  the  resolvent  operator  solution  (73)  can  be 
written  as  the  three-dimensional  convolution  equation 


♦<x>  =  t  e(x)  -  tf"[xT  *0(x)  *  (103) 

The  resolvent  kernel  "drops  off"  with  distance  as  -p  ;  its  effect  is 
thus  localized  to  the  neighborhood  of  the  field  point.  The  resolvent  kernel 
can  thus  be  truncated  at  a  few  wave  lengths  away  from  the  field  point  without 
introducing  much  of  an  error,  and  can  certainly  be  truncated  to  the  complete 
x-domain  P(X)  in  which  a  solution  for  the  field  $(X)  is  sought,  without 
introducing  much  of  an  error  at  all.  Solution  (103)  can  thus  be  represented 
by  the  finite  three-dimensional  convolution 
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(104) 


♦«»  ■  * 6!x) 

PCX) 

This  finite  three-dimensional  convolution  can  clearly  be  evaluated 
numerically  with  the  aid  of  the  fast  Fourier  transform  [22]  algorithm  applied 
to  a  discrete  convolution  [23] ,  for  the  execution  of  which  the  required 
number  of  multiply-add  operations  reduces  to  N  log 2  N  and  the  required  number 
of  storage  allocations  reduces  to  N.  For  the  previously  mentioned  practical 
case  of  N=106,  N  log2  N  ss  2x107,  and  is  thus  executable  on  the  CDC  7600  in 
less  than  five  minutes  [24] . 

For  large  number  of  far-field  data  points  ip ( kc ) ,  it  similarly  becomes 
desirable  to  numerically  evaluate  the  ansatz  (43)  with  the  aid  of  the  fast 
Fourier  transform.  This  is  possible  since  the  ansatz  integral  (43)  is 
essentially  a  Fourier-surface-integral  in  the  observation  transform  k-space; 
i.e.,  since 


k2 


(105) 


where  ds^  is  the  differential  k-space  surface  element,  it  follows  from  (43) 
that 


OCX)  = 


>{i(k)  ds^ 


(106) 


where  the  closed  surface  integration  in  (106)  is  over  a  spherical  k-space 
surface  of  radius  k0.  (The  full  details  of  the  application  of  the  FFT  to 
(106)  will  be  presented  in  subsequent  report (s)  covering  the  follow-on 
contract  to  the  contract  covered  by  this  report) . 


SECTION  X 


FOUR-DIMENSIONAL  FORMULATION  OF  THE  INVERSE  THEORY 


X.l.  INTRODUCTION 


The  objective  of  the  four- dimensional  formulation  of  the  inverse 
problem  is  to  develop  an  exact  inverse  theory  applicable  to  time  dependent 
(i.e.,  moving,  rotating,  and  deforming)  radiators  and  scatterers,  utilizing 
wide-band  short -pulse  Doppler  data.  Some  aspects  of  this  objective,  such  as 
wide-band  inverse  scattering,  can  readily  be  obtained  by  reformulating  the 
frequency  domain  inverse  solutions  into  the  time  domain  with  the  aid  of  the 
time-frequency  Fourier  transform.  However,  other  aspects  of  this  objective, 
such  as  the  utilization  of  Doppler  data,  are  at  best,  extremely  difficult  to 
obtain  by  such  means.  The  source  of  this  difficulty  is  that  all  time 
dependent  solutions  must  be  invariant  under  a  Lorentz  transformation,  which 
is  notoriously  cumbersome  in  three-dimensional  formulations  (particularly 
if  higher  than  velocity  order  effects  such  as  acceleration  and  jerk  are  not 
to  be  neglected)  and  remarkably  simple  in  four-dimensional  formulations. 

A  four-dimensional  relativistically  invariant  formulation  of  the  inverse 
problem  is  thus  a  most  logical  approach. 

A  four-dimensional  scalar  theory,  consistent  with  the  three-dimensional 
scalar  theory  developed  in  the  previous  sections  is  presented  next.  Since 
the  method  of  derivation  of  the  four-dimensional  inverse  integral  equation 
and  its  solution  parallels  closely  the  method  of  derivation  of  the  previously 
presented  three-dimensional  theory,  this  derivation  is  presented  in  a  brief, 
but  complete,  form.  However,  since  again  a  distribution  theory  approach  is 
taken,  this  derivation  is  commenced  with  the  brief  but  concise  rederivation 


of  the  well-known  four-dimensional  formulation  of  Stratton  [25]  of  the  direct 
problem,  using  this  distribution  theory  approach.  The  purpose  of  this 
rederivation  is  not  its  elegance,  but,  as  was  also  the  rationale  for  the 
presentation  of  appendix  I,  its  lending  credence  to  the  distribution  theory 
approach  taken. 


X.2.  FOUR- DIMENSIONAL  FORMULATION  OF  THE  DIRECT  PROBLEM 


Gauss'  theorem  applied  to  the  four-dimensional  vector  field 
(g  o<J>  -  <}>  eg)  yields 


f 


•/. 


da*(g  o<|>  -  <p  Og)  =  I  dr  o*(g  □<}>-$  og) 


(107) 


where  a  is  a  closed  three-dimensional  hyper-surface  in  a  four-dimensional 
hyper- volume  x,  where  the  four-dimensional  scalar  field  <p  and  free-space 
Green's  function  g  satisfy  the  inhomogeneous  wave  equations 


H  =  -  P  (108) 

2g  =  -  6  (109) 


respectively,  and  where 


1 

g  = - -  ;  r  =  x  -  x'  ,  r  =  |r|  .  (110) 

4ir2r2 


Thus,  by  (107)  through  (109), 


^ da*  (g  o<f>  —  4>  og)  =  f  dr  (og*0(j>  +  g  a2$  -  o4>*og  -  <f>  o2g)  (131) 

a  *r 


I  di  (g  d2^  -  <f>  o2g) 

T 


-/ 


dt  g  p  +  :  dx  <j>  6 

T  T 


I 


dr  4>  <$  =  4> ’  ,  for  all  X1  6t 


(112) 


(113) 


(114) 


Thus,  by  (113)  and  (114), 


<j>'  =  I  dt  g  p  + 

da* (g  o<f>  ~  4>  og)  , 

*  T 

J 

a 

(115) 


which  is  the  solution  of  Stratton  [26]  for  the  four-dimensional  direct 
radiation  problem. 


X.3.  FOUR-DIMENSIONAL  FORMULATION  OF  THE  INVERSE  PROBLEM 


Gauss’  theorem  applied  to  the  four-dimensional  vector  field  ( g<J>) 


yields 


da*Q(g<f>)  =  j  dx  o*o(g<j>) 
a  t 

*  j  dr  a 2  ( g4> )  , 

T 


(116) 


(117) 
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which,  by  (108)  and  (109)  yields 


f 


da*(g  o<j>  +  <j>og) 


dx  (g  o<j)  +  <f>  og) 


dr  (Og*a<|>  +  g  o2<{>  +  0<ji‘Dg  + 


dt  (20g*Q<{>  -  g  p  -  <f>  6) 


<j>  o2g) 


(118) 


(119) 


(120) 


If  all  the  non-vanishing  source  distribution  is  in  t  (i.e.,  p(x)=0 
for  all  X  ^  x),  then  by  (115) 


dx  g  p  =  <J>* 


(121) 


Furthermore , , 


I 


dx  4>  6  =  <J>*  ,  for  all  x’  G  t 


(122) 


Thus,  with  the  aid  of  (121)  and  (122),  (120)  yields 


f 


da*(goij>  +  4>og)  =  2  f  dx  og*Oi}>  -  2<f>’ 
a  ■  x 


/■ 


-  J  dx  Og*Di}>  =  -  £  (j>  ( 


<f>'  -  f  dx  og •□$  =  -  i  (|)  da*  (g  o$  +  <f>  og)  , 
x  a 


(123) 


(124) 


which  is  the  four-dimensional  and  integro-differential  equation  representation 
of  the  inverse  problem. 


X.4.  SOLUTION  OF  THE  FOUR-DIMENSIONAL  INVERSE  PROBLEM 


Let  6  be  the  known  four- dimensional  ansatz  to  the  inverse  problem,  i.e., 

6  =  -  i^"da*(gC34>  +  4>  og)  .  (12S) 

a 

The  far-field  ansatz  is  thus  represented  by  a  three-dimensional  hyper¬ 
surface  at  infinity,  and  the  resulting  four-dimensional  domain  of  integration 
over  t  in  (124)  is  thus  over  all  of  four-space.  For  such  a  far-field  ansatz, 
(124)  thus  becomes 


-J 


dx  og*D<j>  =  0' 


(126) 


which,  by  virtue  of  (110) ,  in  four-dimensional  cartesian  coordinates  is 


*(X') 


00 


—CO 


ogCx'-x)'0(j>(x)  d4x  =  e(x') 


(127) 


which  can  be  represented  by  the  four-dimensional  convolution  equation 


$  -  og?  a<}>  =  e 


(128) 


The  four-dimensional  Fourier  transform  of  (12S)  thus  yields  (with  the 
aid  of  the  appropriate  differentiation  rule  □  -*-*■  ik)  in  the  four-dimensional 
transform  k-space 


$  +  k2  G  $  =  0 


(129) 
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where 


$  "4— >■  (j) 

0  ++  0 

> 

(130) 

G  *+  g 

• 

(131) 

Similarly,  (109)  yields  in  four- dimensional  k-space 


(132) 


Combining  (129)  with  (132)  thus  yields 


i  =  f  0 


(133) 


which,  back  in  x-space,  yields 


4>  =  a  © 


(134) 


which,  with  the  aid  of  (125),  yields  the  desired  solution 


(135) 


where  the  domain  of  integration  is  a  closed  three-dimensional  hyper-surface 
at  infinity. 
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It  is  most  remarkable  that  the  resolvent  operator  of  the  four¬ 
dimensional  inverse  integral  equation  (126)  is  pure  local -algebraic  (i.e., 
of  the  form  (134)),  vis-a-vis  the  resolvent  operator  of  the  three-dimensional 
inverse  integral  equation  (25),  which  is  a  global  integral  operator  (i.e., 
of  the  form  (73)).  This  simplicity,  which  occurs  in  four-dimensional  space 
only,  is  intimately  connected  with  another  remarkable  property  unique  to 
four-dimensional  space;  namely,  that  the  functional  form  of  the  Green's 
function  in  k-space  is  invariant  to  the  dimensionality  of  the  space,  which  is 
not  the  case  in  x-space;  and,  that  only  in  four- dimensional  space,  is  the 
functional  form  of  the  Green's  function  the  same  both  in  x-  and  k-space  [27]; 
i.e. , 

G(k)  =  p  ^  “2  =  9{x>  (!36) 

in  four-dimensional  space  only. 

X.5.  CONCLUDING  REMARKS 

The  extension  of  the  preceding  four- dimensional  scalar  field  inverse 
theory  to  vector  and  tensor  fields,  appropriately  applicable  to  electro¬ 
magnetic  fields,  can  be  accomplished  by  means  somewhat  analogous  to  those  of 
sect,  V. 

The  Lorentz- invariant  four- dimensional  reformulation  of  Maxwell's 
equations  for  the  total  charge  and  current  densities  in  terms  of  the 
fundamental  electromagnetic  fields  (see  appendix  II,  eqn.  (32)  -  (35)), 
which  is  needed  for  the  previously  mentioned  extension  of  the  inverse  theory, 
can  be  accomplished  by  most  conventional  means. 

A  four-dimensional  Lorentz-invariant  generalized  boundary  condition  - 
constitutive  equation  (see  appendix  II,  eqn.  (38)  and  (49)),  which  is 


'  5536518555 
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essential  to  the  final  objectives  of  a  four-dimensional  inverse  theory  (see 
sect.  VI)  is  similarly  obtainable. 

The  various  permutations  and  combinations  of  the  special  case 
solutions  (and  their  details)  for  wide-band  short -pulse  data,  monostatic- 
bistatic  data,  Doppler  data,  and  incomplete  ansatz  data,  etc.,  can  be 
extracted  from  the  general  solution  by  an  appropriate  choice  of  the  four¬ 
dimensional  geometry,  the  coordinate  system,  and  the  domain  of  hyper-surface 
integration  implied  by  the  ansatz  (125)  and  the  general  solution  (135);  i.e., 
by  an  appropriate  choice  of  the  "world-line"  of  the  observer  relative  to  the 
radiator  or  scatterer,  or  vice-versa. 

The  work  outlined  in  the  preceding  four  paragraphs  is  currently  being 
undertaken  under  an  appropriate  part  of  a  follow-on  contract  to  the  contract 
covered  by  this  report,  and  will  be  covered  in  detail  in  subsequent  report (s) 
to  this  follow-on  contract. 

Since,  as  was  the  case  with  the  three-dimensional  inverse  solution  (see 
sect.  IX),  the  four- dimensional  inverse  solution  (135)  can  be  formulated  as  a 
convolution;  it  can  thus  similarly  be  evaluated  numerically  with  the  aid  of 
the  fast  Fourier  transform  applied  to  its  equivalent  discrete  convolution 
representation.  However,  since  such  a  convolution  is  four-dimensional,  the 
size  of  practical  data  (e.g.,  < 64 ) 4  =  1.7x107)  will  utilize  to  the  limit  the 
size  and  speed  of  the  biggest  and  fastest  existing  computer  [28] •  Next 
generation  computers,  such  as  Star  or  Illiac,  will  not  be  taxed  to  the  limit 
of  their  size  and  speed  by  four-dimensional  problems  of  practical  magnitude. 

An  alternative  practical  means  of  evaluating  (135)  is  by  a  hard-wired  Fast 
Fourier  Analyze?  (i.e.,  by  special  purpose  computer,  vis-a-vis  the  software 
compiled  FFT  on  a  general  purpose  computer) .  A  further  alternative  might  be 
a  Fourier  analog  optical  processor. 

In  previous  attempts  at  a  formal  definition  of  the  inverse  problem 
(see  sect.  VI),  the  seemingly  logical  classification  of  the  direct  vs,  the 
inverse  problem  in  terms  of  the  unknown  being  the  effect  vs.  the  cause 
respectively,  has  purposely  been  avoided  for  reasons  that  will  be  discussed 
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next.  The  topological  distinction  between  the  three-dimensional  direct  and 
inverse  problem  has  previously  been  shown  to  be  that  of  whether  the  unknown 
fields  and  sources  are  on  different  or  the  same  side  of  the  closed  surface 
ansatz  integrals  (^ds*(g  V<j>  $  4>  Vg)  respectively.  This  topological 
distinction  clearly  carries  over  into  four-dimensions  for  the  closed  hyper¬ 
surface  ansatz  integral  <^da‘(goi}>  ?  <J>og).  However,  because  of  the  temporal 
nature  of  the  added  fourth  dimension  and  the  (lorentz  transformation  invariant) 
causality  principle,  it  is  now  clear  that  no  Lorentz -invariant  distinction 
between  the  direct  and  inverse  problems  can  be  made  on  the  basis  of  cause 
and  effect. 


APPENDIX  I 


SIMPLIFIED  AND  UNIFIED  REDERIVATION  OF  THE 
INTEGRATION  OF  THE  FIELD  EQUATIONS  FOR  THE  DIRECT  PROBLEM 


Gauss*  theorem.,  when  applied  to  the  vector  .field  (<|>  Vg  -  g  V<f>),  yields 


f 


[ 


dS*(<i>  Vg  -  g  V$)  =  I  dv  V*  ( <f>  Vg  -  g  V4>) 


CD 


provided  C<j>  Vg  -  g  V$}  is  continuous  and  twice-differentiable  in  v  and  on  s. 
These  continuity  conditions  can  be  totally  dispensed  with  if  4>  and  g  and, 
hence,  < <f>  Vg  —  g  V<j>),  are  taken  as  distributions  [29];  specifically,  if  <{>  is 
taken  as  a  field  satisfying  the  inhomogeneous  Helmholtz  equation 


V2$(X>  +  k2  $(X)  =  -  p(X) 


(2) 


and  g  is  taken  as  the  associated  Green’s  function 


g(x|x‘) 


ikr 


4nr 


(3) 


r  =  x  -  x* 


(4) 


r  = 


(5) 
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which  satisfies  the  inhomogeneous  differential  equation 


V2g(X)  +  k2  g(X)  =  -  6(X)  , 


and  where  both  9  and  g  satisfy  the  radiation  condition  at  infinity. 


The  right-hand  side  volume  integral  in  (1)  reduces  to 


J  dv  V-  ( 4>  Vg  -  q  V9)  =  J  dv  < Vg  +  V2g  -  Vg •  V<}>  -  g  V2$) 


=  I  dv  (4  V2g  -  g  V2tf>)  , 
**  v 


which  thus  yields  for  Cl) 


<j>dS’(<}>  Vg  -  g  7$)  =  J'  dv  (<j>  V2g  -  g  V2<J>)  , 


which  is  known  as  Green’s  second  identity  in  Green's  theorem. 


With  the  aid  of  (2)  and  (6),  (9)  reduces  to 


^*ds*(9  Vg  -  g  V<J>)  =  J" dv  [9(-k2g-6)  -  g(-k2<J>-p)] 


=  J"  dv  (-k2g9  -  5<f>  +  k2g4>  +  gp) 


-( 


dv  (-5$  +  gp) 


(10) 


(ID 


(12) 


^ ds » ( 4s  Vg  -  g  V<j>) 
s 


dv  6  <{> 


v 


+ 


dv  g  p 


(13) 


If  the  field  point  X'  is  taken  as  inside  the  volume  v  bound  by  the 
surface  s,  then  the  first  volume  integram  on  the  right-hand  side  of  (13) 
yields  by  the  very  definition  of  the  delta-function 


/ 


dv  o(x-x’)  4<x)  =  4 >(xf)  ;  x’  e  v  . 


(14) 


Similarly,  if  the  field  point  X1  is  taken  as  outside  this  volume  v. 


then 


/ 


dv  5(X-XM  <j>(X)  =  0  ;  X'  $  v  . 


(15) 


And  similarly  again,  if  the  field  point  Xf  is  on  the  surface  s 
bounding  the  volume  v,  then 


I 


dv  S(X-X’)  4>(X)  = 


ft(X’ ) 
4ir 


$(X’> 


x>  e  s 


(16) 


where  n(X’)  is  the  internal  solid  angle  subtended  by  s  on  X',  Furthermore, 
if  the  curvature  of  s  at  X'  is  finite  and  non-singular,  then  ft(X*)  =  2n;  then 


f 


dv  SCx-X1 )  4><X)  =  i  <Kxf) 


(17) 


If  the  internal  solid  angle  f»£x’)  subtended  by  s  on  X*  is  consistently 
generalised  for  the  field  point  X*  not  on  s  .e? 
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(18) 


!4tt  ,  X1  Sv 

o  ,  X'  ^  V  ; 


which  is  geometrically  most  reasonable,  (14)  through  (16)  are  representable 
by 


j dv  6(X-X* 


)  *(X)  = 


n(x*  > 


+(X’>  for  all  X*. 


With  the  aid  of  this  (19) ,  (13)  yields 


*'  8  dv  g  p  +  (j)ds*(g  V(J>  -  $  Vg)  , 


which  is  a  unified  representation  of  the  Kirchhoff  integration  of  the  field 
equations,  valid  for  all  field  points  (i.e.,  field  points  inside,  on,  and 
outside  s) . 

A  physical  interpretation,  as  well  as  a  rigorous  proof  of  the  physical 
meaning  of  both  the  volume  and  the  surface  integrals  in  (20)  are  of  interest 
next.  To  this  end,  consider  a  closed  surface  sj  totally  inside  another 
closed  surface  s2  (see  fig.  12),  a  source  distribution  pj(X)  totally  localized 
in  the  volume  vj  bound  by  the  closed  surface  si  (i.e.,  pj(X)  =  0  for  all  X 
outside  si),  and  another  source  distribution  p2  (distinct  from  p*)  totally 
localized  in  the  intervening  volume  v2  band  by  the  closed  surfaces  si  and  s2 
(i.e.,  p2(X)  =  0  for  all  x  inside  si  and  outside  s2).  Furthermore,  let  the 
field  point  x'  be  outside  the  closed  surface  and  inside  the  closed  surface 
s2j  i.e,,  in  the  volume  v2.  It  thus  follows  that  if  the  closed  surface  s 
in  (20)  is  taken  as  the  closed  surface  S!  and  s2,  then  by  (18),  (20)  yields 


-  I  dv  g  p2  + 


£ ds*  (g  7$  - 


*  Vg)  +  dS* (g  V<{>  -  <t>  Vg) 
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(23) 


i’  =  J  dv  g  Pj 


or  the  closed  surface  integral 


=  ^*dS‘(g  V^  -  <\>1  Vg) 


(24) 


That  these  volume  and  surface  integrals  (23)  and  (24)  are  indeed 
equivalent  can  also  be  shown  by  taking  the  closed  surface  s  in  (20)  as  the 
closed  surface  sj  only;  then,  by  (18),  (20)  yields 


■/ 


dv  g  Pi 
vi 


-  jj>  ds*  (g  V<j>j  -  Vg) 
si 


Q.E.D. ; 


(25) 


(where  the  (-)  sign  in  front  of  the  surface  integral  in  (25)  is  du3  to  the 

* 

fact  that  by  Gauss’  theorem  the  surface  unit  vector  s  must  always  be  chosen 
as  outwards  from  the  domain  of  volume  integration;  see  fig.  12). 

The  physical  meanings  and  interpretations  of  (23)  and  (24)  respectively 
are  clearly  that  of  the  sm  (volume  integral)  of  the  action-at-a-distance- 
field  due  to  a  spatial  (volume)  source  distribution,  and  its  closed  surface 
integral  equivalent ,  provided  the  field  point  X’  and  the  source  distribution 
Pl  are  outside  and  inside  of  the  closed  surface  si  respectively. 

Next,  consider  the  field  point  X"  chosen  as  inside  the  closed  surface 
Si  (see  fig.  12);  thus,  if  the  closed  surface  s  in  (20)  is  tasen  as  the 
closed  surface  si,  then  by  (18),  (20)  yields  for  the  field  <}>”  at  this  field 
point  X" 


!>"  =  J*dv  g  p!  -  (j>( 


<j>"  =  3  dv  g  pi  -  (J)  ds •  ( g  v$  -  $  vg) 

vi  -  Sl 
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(26) 


(where  the  (-)  sign  in  front  of  the  surface  integral  in  (26)  is  again  for 
the  same  reasons  as  stated  subsequent  to  (25)). 

If,  however,  the  closed  surface  s  in  (20)  is  taken  as  the  closed 
surface  si  and  S2,  then  by  (18),  (20)  yields  (for  the  field  point  still 
chosen  as  X") 


0  =  f  dv  g  p2  +  £ 
*v2 


ds*(g  '7<|>  -  <j>  Vg) 
si 


(27) 


It  thus  follows  from  the  previously  developed  physical  meaning  and 
interpretation  of  the  volume  integral  over  a  source  distribution,  and 
comparison  of  (26)  and  (27),  that  the  physical  meaning  and  interpretation  of 
the  closed  surface  integral  in  (26)  is  that  of  the  (negative  of  the) 
equivalent  of  the  action-ai-a-distance- field  due  to  a  spatial  source 
distribution ,  provided  the  field  point  X"  and  the  source  distribution  p2  are 
inside  and  outside  of  the  closed  surface  Si  respectively  (the  (-)  sign  again 
having  arisen  for  the  previously  stated  reasons) . 

Since  in  the  previous  arguments  the  closed  surface,  the  field  point, 
and  the  nonvanishing  source  distribution,  were  chosen  topologically  com¬ 
pletely  arbitrarily,  with  the  sole  exception  that  the  field  point  and  the  non¬ 
vanishing  source  distribution  are  to  be  on  opposite  sides  of  the  closed 
surface,  it  follows  that  if  the  surface  unit  vector  is  consistently  chosen 
as  pointing  from  the  field  point  to  the  source  distribution,  then  the 
physical  meaning  and  interpretation  of  the  surface  integral  (^jds*(g  V<}>  -  $  Vg) 
is  ahsxzys  that  of  the  fields  produced  at  a  field  point  on  ONE  SIDE  of  the 
surface  by  the  total  source  distribution  located  on  the  OTHER  SIDE  of  the 
surface,  and  totally  equivalent  to  the  volume  integral  J  g  p  dv  over  this 
same  source  distribution. 

This  property  of  relating  fields  and  sources  on  opposite  sides  of  the 
Kirchhoff  surface  integral  is  its  fundamental  topological  property,  when 
applied  to  direct  field  problems j  it  is,  however,  totally  useless  when 
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applied  to  inverse  field  problems;  since  the  latter  problem  is  topologically 
characterized  by  unknown  fields  and  sources  on  the  same  side  of  a  closed 
surface.  The  surface  integral  of  interest  to  the  inverse  field  problem  is 
thus  a  surface  integral  relating  fields  and  sources  on  the  same  side  of  the 
surface;  relating  them  at  least  in  such  a  fashion  as  to  render  the  surface 
integral  a  useful  ansatz  to  the  inverse  problem.  As  shown  in  sect,  II,  such 
a  surface  integral  is  (^>ds*(g  +  4>  Vg). 

Replacing  the  scalar  field  Gauss'  theorem  (1)  with  its  vector 
equivalent 


f 


J 


dS*(VgxF  -  g  V*F)  =  I  dv  Vx(Vgxp  ~  9  VxF) 


(28) 


and  applying  it  analogously  twice  to  the  electromagnetic  vector  field  and 
wave  equations  (see  appendix  II),  vis-a-vis  the  scalar  Helmholtz  field 
equation  (2),  where  the  vector  field  F  represents  the  electric  and  magnetic 
vector  fields  respectively,  similarly  yields  the  familiar  Stratton  and  Chu 
[30]  or  Franz  [31]  integral  representation  of  Maxwell's  equations. 

As  in  the  case  of  the  scalar  field,  the  surface  integral  in  (28) 
relates  fields  and  current  density  distributions  (the  source  distributions) 
of  opposite  sides  of  this  surface;  whereas  the  surface  integral  <^>dSxvx(gF) 
relates  fields  and  current  densities  which  are  on  the  same  side  of  the 
surface. 


APPENDIX  U 


A  REFORMULATION  OF  THE  ELECTROMAGNETIC  FIELD  EQUATIONS. 


II. 1.  THE  CONVENTIONAL  FREE  CHARGE  AND  CURRENT  DENSITIES  REPRESENTATION 


The  conventional  representation  of  Maxwell's  equations  is 


7-D  =  „f  (1) 

V*B  =  0  (2) 

-  -  f 

7xH  =  0f+f  ,  (4) 


where  p j.  and  J j.  are  the  free  charge  and  current  densities  respectively  (also 
referred  to  by  some  authors  as  the  true  or  unbound  charge  and  current  densities) ; 
E  and  B  are  the  electric  and  induction  fields  respectively,  which  are  funda¬ 
mental  fore©  fields;  and  B  and  H  are  the  displacement  and  magnetic  fields 
respectively,  which  axe  artificially  derived  fields.  These  artificially 
derived  fields  are  related  to  the  fundamental  fields  by 


D  =  Go 

E  +  Pi 

(5) 

H  =  1 

Vo 

B  -  M  , 

(6) 

where  and  M  are  the  polarization  (of  the  bound  charges)  and  magnetization 
fields  respectively.  These  polarization  and  magnetization  fields  are  the 
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contributions  to  the  fundamental  fields  produced  by  local  charge  and  current 
densities  in  a  material  medium.  e<>  and  y0  are  the  permittivity  and  perme¬ 
ability  of  free  space  respectively,  collectively  referred  to  as  the  inductive 
capacities  of  free  space. 

In  linear  isotropic  media,  these  polarization  and  magnetization  fields 
are  relatable  to  the  electric  and  magnetic  fields  respectively  by  the 
constitutive  equations 

h  =  Xe  E  C7) 

M  =  xm  H  ,  '  (8) 

where  xg  and  are  the  electric  and  magnetic  susceptibilities  of  the  linear 
isotropic  medium  respectively,  collectively  referred  to  as  the  susceptibilities 
of  the  medium. 

For  such  linear  isotropic  media,  with  the  aid  of  (5)  through  (8),  the 
relationship  between  the  derived  and  fundamental  fields  yields  the 
conventional  constitutive  equations 

D  =  e  E  (9) 

B  =  u  H  5  (10) 


where 


E  =  K  Eo 

e 


K  =  1  +  X  t 

e  Ke  * 


I' 


u  . 


=  1  +  *m  ' 


where  e  and  y  are  the  permittivity  and  permeability  of  the  medium  respec¬ 
tively,  collectively  referred  to  as  the  inductive  capacities  of  the  medium; 
and  Kg  and  are  the  relative  permittivity  (or  dielectric  constant)  and 
relative  permeability  of  the  medium  respectively,  collectively  referred  to 
as  the  specific  inductible  capacities  of  the  medium. 

For  linear  isotropic  conducting  media,  Ohm's  law  yields  the 
additional  constitutive  equation 


where  a 


f 


is  the  conductivity  (of  free  charges  of  the  medium. 


(15) 


For  linear  isotropic  media  which  are  also  spatially  and  temporally 
homogeneous  (i.e.,  media  for  which  the  inductive  capacities  are  constants 
neither  a  function  of  space  nor  a  function  of  time),  the  conventional 
formulation  of  the  wave  equations  follows  from  (3),  (4),  (9),  and  (10);  i.e.. 


VxvxE  +  ye  -0  =  -  y  |^f  (16) 

vxvxH  +  ye  0  =  Vxj^  _  (17) 


For  the  general  case  of  inhomogeneous  media  (i.e.,  media  for  which 
either  or  both  of  the  inductive  capacities  are  not  a  constant;  namely,  are  a 
function  of  time  and/or  space) ,  the  formulation  of  such  a  wave  equation  is  in 
general  not  possible. 

Ihe  preceding  equations  can  be  classified  into  one  of  three  classes: 

MAXWELL'S  EQUATIONS:  The  first  order  differential  equations  (1) 
through  (4)  relating  the  electric  and  magnetic  fields  to  the  charge  and 
current  densities.  In  these  equations  the  electric  and  magnetic  fields  are 
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intertwined  and  not  separated.  These  equations  are  the  fundamental  equations 
derivable  from  the  fundamental  action-at-a-distance  laws  of  Coulomb,  Faraday, 
and  Biot-Savart. 

THE  WAVE  EQUATIONS:  The  second  order  differential  equations  (16), 

(17)  relating  the  electric  or  magnetic  fields  separately  to  the  current 
densities..  The  essential  and  distinguishing  (from  Maxwell’s  equations) 
property  of  these  equations  is  that  tht  vlectric  and  magnetic  fields  are 
separated,  i.e.,  that  each  wave  equation  contains  one  field  (electric  or 
magnetic)  only.  The  term  wave  equation  is  thus  possibly  a  misleading 
misnomer;  a  more  characteristically  descriptive  name  might  thus  be  the 
separated  equations.  It  should  again  be  noted  that  for  the  general  case 
of  inhomogeneous  media,  the  separation  of  the  electric  and  magnetic  fields 
into  a  set  of  wave  equations  of  the  form  (16)  and  (17)  cannot  be  achieved. 

Maxwell's  equations  and  the  wave  equations,  namely,  all  those 
equations  derivable  from  the  fundamental  action  at  a  distance  laws,  are 
collectively  referred  to  as  the  field  equations. 

THE  CONSTITUTIVE  EQUATIONS:  The  local  algebraic  equations  (7),  (8), 
(9),  (10),  and  (15)  which  relate  the  fundamental  fields  to  the  artificially 
derived  fields,  or  the  current  densities  to  the  fields.  These  equations  are 
a  description  of  the  local  properties  of  the  material  medium,  such  as  Ohm's 
law,  and  are  essentially  non-electromagnetic  in  nature  (namely,  they  do  not 
contain  the  basic  electromagnetic  action-at-a-distance  laws);  in  most  cases, 
these  equations  involve  Newton's  second  law  of  motion  and  the  microscopic 
dynamics  of  the  material  medium. 

It  takes  a  (set  of)  field  equations  and  a  (set  of)  constitutive 
equations  to  completely  specify  any  electromagnetic  problem.  (The  consti¬ 
tutive  equations  often  take  the  form  of  a  boundary  condition  to  the 
differential  field  equations;  the  boundary  conditions  determining  the 
specific  unique  solution  to  the  differential  equations) . 


100 


II. 2.  THE  DIFFICULTIES  WITH  THE  CONVENTIONAL  FORMULATION 


The  difficulty  with  the  conventional  formulation  of  the  field 
equations  is  essentially  the  severely  restricted  applicability  of  the 
resulting  wave  equation  to  the  interior,  and  not  the  surface,  of  homogeneous 
media  only. 

It  thus  immediately  follows  that  the  conventional  formulation  of  the 
integral  representation  of  the  field  equations  is  also  valid  for  the  interior, 
and  not  the  surface,  of  homogeneous  media  only;  i.e.,  the  complete  domain 
of  volume  and  surface  integration  must  be  in  the  interior  of  the  homogeneous 
medium. 


The  seriousness  of  this  limitation  can  best  be  illustrated  by  the 
case  of  a  finite  sized  homogeneous  medium  of  constant  inductive  capacities 
different  from  those  of  empty  space,  imbedded  in  infinite  free  space.  For 
such  a  case,  one  clearly  obtains  two  different  wave  equations  (of  the  form 
of  (16)  and  (17) ,  one  valid  inside  the  medium  and  one  valid  outside  of  the 
medium,  neither  of  which  is  valid  on  the  surface  between  the  medium  and  the 
free  space.  Furthermore,  there  exists  no  single  wave  equation  valid  for  all 
space.  The  conventional  formulation  of  the  integral  representation  of  the 
field  equations  is  thus  valid  only  if  the  complete  domain  of  volume  and 
surface  integration  is  completely  in  the  interior  of  such  a  medium;  this 
integral  representation  is  thus  not  valid  for  such  problems  as  the  external 
scattering  by  such  media. 

For  the  general  case  of  arbitrarily  inhomogeneous  media,  a  wave 
equation,  though  much  more  complicated  than  the  form  (16)  and  (17)  is 
formable;  however,  no  appropriate  Green's  function  is  known  for  such  a  wave 
equation.  Without  such  a  Green's  function,  the  integration  of  the  field 
equations  into  an  integral  representation  is  thus  not  possible. 


mmmm. 
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II. 3.  THE  TOTAL  CHARGE  AND  CURRENT  DENSITIES  REPRESENTATION 


The  source  of  the  difficulties  with  the  conventional  formulation  of 
the  field  equations  lies  in  their  provincial  view  of  matter,  charge,  and 
current,  and,  hence,  the  fields.  Maxwell’s  equations  contain  only  the  free 
charge  and  current  densities,  hence,  the  artificially  derived  fields  in 
addition  to  the  fundamental  fields.  The  presence  of  the  artificially  derived 
fields  in  Maxwell’s  equations  furthermore  injects  into  the  field  equations 
some  of  the  properties  of  the  medium,  which  rightfully  should  be  in  a 
separate  set  of  constitutive  equations.  The  distinction  between  free  and 
bound  charge  and  current  densities  is  artificial  and  a  source  of  considerable 
difficulty  when  dealing  with  non-classical*  non-linear  media  such  as  plasmas 
and  semiconductors.  Furthermore,  all  accelerating  charges,  whether  free  or 
bound,  radiate  electromagnetic  fields.  It  thus  stands  to  reason  that  a  more 
fundamental  formulation  of  Maxwell’s  equations  would  contain  only  the  total 
charge  and  current  densities  and  only  the  fundamental  fields,  and  thus  be  free 
of  the  artificially  derived  fields  and  the  artificial  distinction  between  free 
and  bound  charge  and  current  densities.  The  properties  of  the  medium  implicit 
in  the  artificially  derived  fields  should  then  only  appear  in  a  separate  set 
of  constitutive  equations.  Such  reasoning  is  further  enhanced  by  the  point 
of  view  of  an  observer  of  electromagnetic  radiation,  to  whom  in  principle  no 
distinction  exists  between  electromagnetic  radiation  emanating  from  accelera¬ 
ting  free  or  bound  charges,  and  who  measures  such  radiated  fundamental  fields 
at  an  observation  point  with  the  aid  of  the  properties  of  the  medium  in  which 
the  observation  point  is  imbedded. 

The  sought  reformulation  of  Maxwell’s  equations  is  thus  a  formulation 
invariant  to  the  specific  properties  of  the  medium  and  the  mechanistic  nature 
of  the  free  and  bound  charges  in  the  medium;  namely,  a  formulation  stating 
the  universal  laws  of  the  electromagnetic  theory  only.  All  the  specific 
properties  of  the  medium  would  thus  be  relegated  solely  to  a  set  of 
constitutive  equations . 
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Such  an  invariant  formulation  of  Maxwell's  equations  can  be  derived 
directly  from  the  basic  aotion-at-a-distanoe  laws  of  Coulomb,  Faraday,  and 
Biot-Savart,  and  an  argument  analogous  to  Maxwell's  argument  that  yielded 
the  displacement  current,  utilizing  an  equation  of  continuity  for  the  conser¬ 
vation  of  the  total  charges  only.  However,  for  the  sake  of  showing  consistency 
with  the  conventional  formulation  of  Maxwell's  equations,  such  an  invariant 
formulation  will  be  developed  next  directly  from  the  conventional  formulation 
and  the  definitions  of  the  artificially  derived  fields. 

Eliminating  the  artificially  derived  fields  from  (1)  through  (4)  with 


the  aid  of  (5)  and  (6)  yields 

VE  =  i  to,  - 

V.Pj,) 

(18) 

v*B  =  0 

(19) 

VxE  =  -  & 

V  11  9t 

(20) 

VxB  =  Uo  (0^  + 

+  vxM)  +  yoCo  ^ 

(21) 

From  the  very  definitions  of  the  polarization  and  magnetization  fields 
associated  with  the  material  media,  it  follows  that 


n  -  -  ”f 

(22) 

— 2>  =  J 

9t  p 

(23) 

VxM  =  J  , 

m 

(24) 

where  is  the  bound  polarization  charge  density;  namely,  the  net  bound  charge 
density  produced  by  the  polarization  field  in  a  material  medium;  and  and 
are  the  bound  polarization  and  magnetization  current  densities  respectively; 
namely,  the  net  bound  current  densities  produced  by  the  polarization  and 
magnetization  fields  respectively  in  a  material  medium. 
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Since  the  total  charge  and  current  densities  p  and  d  respectively  is 
the  sun  of  all  the  free  and  bound  charge  and  current  densities  respectively, 
irrespective  of  the  mechanistic  nature  of  these  bound  charge  and  current 
densities,  it  follows  that 


p  =  p /  "  7*Pi 


d  =  <3y»  +  +  vxM 


The  invariant  formulation  of  Maxwell's  equations  in  terms  of  the  total 
charge  and  current  densities  thus  becomes 


v*E  =  -  p 

£o 


v*B  =  0 


vxE  =  - 


VxB  =  y0  d  +  UoEp  tz 

du 


It  now  becomes  convenient  to  redefine  the  magnetic  field  H  as  a  funda¬ 
mental  field,  related  to  the  fundamental  induction  field  B  for  all  media  as 


H  =  -  B  . 

Uo 


Maxwell's  equations  (27)  through  (30),  in  terms  of  this  redefined 
fundamental  magnetic  field  thus  become 


v-E 

1 

=  —  p 

So 

v*H 

=  0 

vx£ 

3H 

’  ~  y°  9t 

vxH 

_  ,  ,  dE 

-  0  +  e«  3t  • 

t 
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The  invariant  formulation  of  the  wave  equations,  no  longer  restricted 
to  linear,  isotropic,  and  homogeneous  media,  but  universally  valid  for  all 
media,  now  becomes 

(?6) 

(37) 

The  three  constitutive  equations  (7) ,  (8) ,  and  (15) ,  with  the  aid  of 
C23)  and  (24),  thus  immediately  yield  the  single  constitutive  equation  for 
the  total  current  density 

0  =  Oj.  E  +  so  E)  +  H)  .  (38) 

The  boundary  condition  for  the  surface  of  perfect  conductors,  i.e., 

K  =  n*H  ,  (39) 

where  K  is  the  surface  current  density  and  n  is  the  outward  surface  unit 
vector,  presents  a  special  problem  that  can  be  resolved  by  introducing  the 
magnetic  vector  surface  conductivity  distribution  density  (dv}  such  that, 
since  0  dv  =  K  ds,  the  boundary  condition  (39)  becomes  the  constitutive 
(distribution)  equation 

J  =  |xH  .  (40) 

With  the  aid  of  (38),  it  is  thus  possible  to  introduce  the  generalized 
constitutive  equation 


* 

■5 

i 


The  equation  of  continuity  for  the  conservation  of  the  total  charges 

thus  is 

7-0  +  |f  =  0  .  (42) 

9t 

Hie  time-reduced  form  of  Maxwell's  equations  (32)  through  (35)  for 
monochromatic  fields  in  terms  of  the  total  charge  and  current  density  thus  is 


v*E  = 

1 

—  p 

Co 

7*K  * 

0 

u 

UJ 

X 

o 

-  itiiyo  H 

7xH  = 

d  +  iue©  E  . 

With  the  aid  of  (44) ,  the  time-reduced  form  of  the  wave  equations  (36) 
and  (37)  for  monochromatic  fields  in  terms  of  the  total  current  density  thus  is 


VxVxE  -  k2  E  =  -  iwy©  d 
V2H  +  k2  H  =  -  VxJ 


The  time-reduced  form  of  the  generalized  constitutive  equation  (41) 
for  monochromatic  fields  for  the  total  current  density  thus  is 


d  =  (cy  +  iue©  Xg)E  +  \  V  +  i 


And  the  time-reduced  form  of  the  equation  of  continuity  for  the 
conservation  of  total  charges  thus  is 

V*d  =  -  in  p  . 

A  universally  valid  scalar,  vector,  and  Hertz  potential  theory  for 
the  fundamental  fields  only,  in  terras  of  the  total  charge  and  current 
densities,  can  similarly  and  consistently  be  developed. 
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13.  ABSTRACT 


PART  I:  A  three-diemsnional  electromagnetic  Inverse  Scattering  Identity,  based  on  the  Physical 
Optics  approximation,  is  developed  for  the  monostatic  scattered  far  field  cross  section  of  perfect  con¬ 
ductors.  Uniqueness  of  this  inverse  scattering  identity  is  proven.  This  identity  requires  complete 
scattering  information  for  all  frequencies  and  aspect  angles.  A  non-singular  integral  equation  is  dev¬ 
eloped  for  the  arbitrary  case  of  incomplete  frequency  and/or  aspect  angle  scattering  information.  A 
general  closed  form  solution  to  this  integral  equation  is  developed,  which  yields  the  shape  of  a  scatterer 
from  such  incomplete  information.  A  specific  practical  radar  solution  is  presented.  The  resolution  of 
this  solution  is  developed;  yielding  short-pulse  target  resolution  radar  system  parameter  equations. 

Results  of  the  three-dimensional  numerical  reconstruction  of  a  sphere  and  a  cylinder  from  a  variety  of 
aspect  angle  and  frequency  band  limited  cross  section  data  are  presented.  The  merits  of  this  solution 
over  the  conventional  synthetic  aperture  radar  imaging  technique  are  discussed. 

PART  II:  The  general  inverse  scattering  and  radiation  problem  associated  with  the  three-dimensional 
inhomogeneous  scalar  field  Helmholtz  wave  equation  is  formulated  as  a  Fredholm  integro-differential  equation 
of  the  second  kind.  The  far-field  inverse  integro-differential  equation  is  solved  in  closed  form  with  the 
aid  of  a  single  resolvent  integral  operator,  which  can  be  readily  evaluated  numerically  with  the  aid  of  the 
Fast  Fourier  transform  algorithm.  The  inverse  integro-differential  equation  and  its  solution  are  then 
generalized  to  the  reduced  vector  wave  equation  resulting  from  Maxwell's  equations.  A  formal  statement  of 
the  inverse  problem  is  presented.  It  is  shown  that  the  first  order  Neumann  series  solution  of  the  inverse 
integro-differential  equation  as  well  as  the  fiTst  order  term  of  its  exact  solution  represent  the  physical 
optics  approximation  and  the  equations  governing  synthetic  microwave  holography.  The  (analogous)  four- 
dimensional  inverse  integro-differential  equation  and  its  closed  form  solution,  applicable  to  Doppler-con¬ 
taminated  fields  radiated  by  moving  scatterers,  is  developed. 
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